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ABSTRACT 



The two principal existing methods of calculating axially- 
symmetric compressible flov/ in turborr.achincs are: (1) a simplified 
one -dimensional analysis, and (2) numerical methods using the com- 
plete or linearized flow equations. The first is not satisfactory for 
multi-stage turbines with appreciable wall divergence; the second is 
very tedious and time consuming. The purpose of this investigation 
is to extend the approximate methods, successfully used in calcu- 
lating incompressible flow in compressors with constant blade height, 
to the analysis of compressible flow in turbomachines with variable 
blade height. Assuming that the blades can be completely defined by 
the exit flow angle, and neglecting the influence of downstream blades, 
the analysis is made considering the flow between successive blade 
rows only. Ivith these restrictions, subsonic and isentroplc super- 
sonic flow patterns can be determined for arbitrary boundary shapes 
as long as separation does not occur. Average losses can be ac- 
counted for by the use of a polytropic law, and the effect of radial 
variations in stagnation temperature can be included without diffi- 
culty. Examples illustrating the flexibility and practical value of the 
iteration method, and the rapid convergence of successive solutions 



are given. 
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I. INTRODUCTION 

In the modern aircraft gas turbine, the flow is necessarily 
in the transonic range for two reasons: (1) to keep the weight and 
size of tlie machine at a minimum, and (2) to fit the mass flow and 
pressure ratio of tlie compressor. Earlier turbojet engines had single 
stage transonic turbines, but present engines employ multistage tur- 
bines with rapid channel divergence in order to use effectively the 
increased pressure ratios and keep the axial velocities compatible 
with the blade speed. In these types of machines, the flow at the 
mean blade height is mainly axial, and the aspect ratio of the blades 
is usually above unity for weight saving considerations. As long as 
the channel wall divergence is not too great so that separation can 
be avoided, there is no reason for low efficiency, particularly in 
view of the favorable pressure gradient through the turbine. 

Two distinct design problems occur in practice which Marble 
(Reference 15) has classified as the direct and inverse problems of 
turbomachine design in analogy with the corresponding classical 
problems in the theory of finite wings. The direct, or "off-design” 
problem occurs when tiie blade speed, blade shape, boundary configu- 
ration, and fluid state ahead of all blades is given, and the three- 
dimensional velocity field, blade loading, and distribution of energy 
are to be determined. The inverse problem occurs In the original 
design of turbomachinery, when the blade loading, blade speed, 
boundary configuration, and fluid state ahead of all blades is pre- 
scribed, and the three-dimensional velocity field, blade shape, and 
distribution of energy in the field are to be determined. 

The exact analysis of viscous, compressible flow through an 



a:naL turbomachine will probably ncvei* be poooible because of the 
non-linear character of the equations governing the flow. Hov/cvcr, 
the design engineer has a great need for an analysis suitable for 
engineering purposes, and may have no use for a lengthy, compli- 
cated and costly procedure even if it were available. He needs an 
extremely flexible metliod of analysis which, while being as simpli- 
fied as practicable, still produces information accux*atc enough for 
effective design woi*k. 

The equations governing three-dimensional, compressible 
channel flow are greatly simplified i£ terms involving viscosity arc 
neglected. In viev/ of the fact that viscous effects actually are very 
slight except near the boundaries, and as long as separation is avoided, 
the assumption of zero viscosity seems justified. Viscosity effects 
could be accounted for as an added correction to tlie basic solution. 
However, with the present day turbines v/hich operate in the high 
subsonic range, it is obvious that compressibility effects can not 
be neglected. 

One of tlie principal methods available to the design engineer 
for the analysis of axially-symmetric, compressible flow in turbo- 
machines is a simplified one -dimensional analysis. This method is 
used extensively in preliminary design considerations and is simple 
and straightforward. Since only the exit flow angles from the blade 
rows at the mean blade height are considered, the radial distributions 
of velocity and energy are not determined. Although this analysis 
gives reasonably accurate results for a single stage turbine with a 
fairly high hub ratio and small channel divergence, the one -dimensional 
method is inadequate for large channel divergence where the radial 
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distribution of velocities becomes significant. It is particularly in- 
adequate for the design of multistage turbines v/ith channel wall di- 
vergence since it is difficult to get the proper balance of worL on the 
individual stages. This work balance is strongly affected by the radial 
distribution of the axial and radial velocities. .?ince the one-dirr.en- 
sional analysis is done at the mean blade height, it gives no informa- 
tion as to the radial distribution of velocities. 

The second principal method of analysis of axially- symmetric 
compressible flow in turbomachines is by mvmerical methods using the 
complete or linearized flow eqiiations to solve for the complete velocity 
and energy field through tlie machine. Since the determination of the 
entire flow field Involves six unknowns at each point in the flow (the 
three components of velocity, the pressure, the temperature, and the 
density), the analysis using these numerical methods is necessarily 
lengthy and extremely complicated and cannot be considered suitable 
from the point of view of the design engineer. Monroe (Reference 7) 
formulated the idealized problem in terms of a stream function for the 
velocities in the meridional plane, and solved the resulting non-linear 
differential equation by a simultaneous application of the relaxation 
technique of Southwell and an iteration process. Wu (Reference 9) has 
utilized various relaxation methods, matrix systems, and finite differ- 
ence schemes to solve the complete hydrodynamic al equations for the 
entire flow field in many types of problems. Vavra (Reference 8), 
Reissner (Reference 10), Goldstein (Reference 11), Wislicenus (Ref- 
erence 13), and many others have employed various techniques to solve 
the complete, or nearly complete, equations by some numerical method. 
Although these numerical methods are not suitable for design purposes. 
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they are quite general and can be used for centrifugal, mixed flow or 
axial flow machines v/ith blades of any aspect ratio. 

In an effort to simplify the problem of solving for the complete 
flow field, a linearized analysis has been provided by Marble (Refer- 
ences 5 and 15) and Fabri (Reference 6). Marble (Reference 15) pro- 
posed a linearized analysis of the problem, satisfactory for high aspect 
ratio blades, to allow an approximate treatment of the general blade row 
with prescribed wall geometry. He assumed that the trailing vorticity 
could be considered to be transported downstream by the mean axial 
flow and not by the perturbation velocities, that is, disturbances in the 
radial and axial velocities are small in comparison with the mean axial 
velocity, and that the blade row is made up of an infinite number of 
blades of finite chord. The analysis for the inverse problem was carried 
out in detail and v/as found to allow a general and reasonably simple so- 
lution. In Reference 5, Marble extended the analysis to the direct prob- 
lem of off-design operation, and to tlic study of mutual interference of 
neighboring blade rows in a multistage axial turbomachine. He provided 
examples of axial and conical flow, and stated that the simple linearized 
solution was sufficiently accurate if the vorticity effects were not large. 
A second order linearization was given to handle problems with greater 
vorticity effects. However, the computation required to find the three- 
dimensional flow in any particular case is rather lengthy using the lin- 
earized solutions, and involves some type of numerical integration. A 
simple exponential approximation was then introduced to simplify the 
computation, and the results using this approximation compared favor- 
ably with the detailed results of the linearized solution. 

The analysis of flow in turbomachinery where tlie channel bound- 
aries are purely cylindrical, or vary only slightly from a cylindrical 
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shape, has been extensively investigated by Raanie (References 1 and 14), 
Traupel (References 2 and 3), Sinnette (References 4 and IS), and others. 

In these analyses, the radial momentum equation used in the treatment is 
of the form, 

/ 6>;o _ 

e dr r (1), 

where the radial force and other momentum terms have been neglected. 

The distribution of axial velocity is approximated by assuming axial sym- 
metry, and since the centrifugal force within the rotating fluid body is bal- 
anced only by the radial pressure gradient, the resulting flow calculated 
in this manner is in reality that which must exist far downstream of the 
blades where radial velocities and accelerations have vanished. No in- 
formation is provided by these analyses on how rapidly the cliange in ve- 
locity' pattern takes place as the fluid passes through the blade row. 

The purpose of tliis thesis is to present an approximate method 
for the analysis of compressible flow in modern axial turbomachines which 
gives sufficiently accurate information for design purposes and yet retains 
the simplicity necessary to be of practical value to the design engineer. 

This is the logical step forward from the one -dimensional method to the 
analysis of "compressible flow in axial turbomachines. Methods which 
have been successfully employed in the incompressible case (Reference 1), 
liave been modified and extended to tlie analysis of compressible flov/ in 
arbitrarily shaped axially-symmetric channels. It is obvious tliat little 
simplification over the so-called "numerical" methods can be obtained if 
the complete flow and energy fields are to be determined in the analysis. 

In this thesis, the flow is analysed at tlie station immediately downstream 
of the blade row by a consideration of tlic upstream conditions, channel con- 
figuration, and blade shapes. The blade trailing edge angle only is consid- 
ered, and the complete flow field through the blade rov/ is not determined. 



The method of analysis presented is extremely flexible and has 
been used to solve the direct and inverse problems, subsonic or super- 
sonic, for arbitrary asiial boundary configurations. In the solution of 
the direct problem, the choice of blade speed, blade spouting angle, 
boundary configuration, and upstream conditions can be arbitrary, and 
the three-dimensional velocity field, blade loading, and distribution of 
energy in the field are determined. In the inverse problem, for a given 
blade loading, blade speed, boxmdary configuration, and fluid state ahead 
of all blades, tlae blade spouting angle, velocity field and energy distri- 
bution can be determined. By a simple change of sign in the equations 
used, the method can be applied equally well to an axial compressor as 
to an axial turbine. Although tliis method is not as general as those used 
by "Wu (Reference 9), for example, it is more flexible, and practical for 
blading of aspect ratio greater than unity, and is sufficiently accurate 
in view of the assumptions used in the treatment. 

The analysis is carried out on the basis of steady, axially- 
symmetric flow with infinitely many blades. The fluid is assumed to 
be compressible, but all other real fluid effects are neglected except 
as they may be approximated by a "polytropic efficiency" and the use 
of the polytropic exponent "n". Ae a result of this assumption, the jjos- 
sibility of separation of flow at the channel walls or along the blades is 
excluded. Boundary layer effects are also not considered in the anal- 
ysis. In the supersonic solution, the occurrence of shock waves is 
excluded. Radial or nearly radial blades are presupposed so that the 
radial force can be neglected, and the radial momentum equation used 
in the treatment has the form. 
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This asauinption seems justified in view of the investigations of -,'^.arlsson 
(Reference 19)» who, by a tv/o -dimensional analysis of incompressible 
flow v/itli an infinite number of blades, showed that the effect of the radial 
force resulting from normal blade twist is of negligible magnitude. 

The effect of downstream blades on conditions upstream is ne- 
glected, and the blade shape is specified by the trailing edge angle only. 
Simply stated, a blade rov/ is assumed whose only effect on the flow is 
to turn it through a specified angle distribution which is a function of 
radius across the channel. An assumption similar to this has been used 
successfully in the analysis of flow through a compressor (References 
1 and 14), and should be even more applicable to turbines where the blade 
solidity is higli. The equations arc particularly suited to the analysis 
of flows with a prescribed radial total entlaalpy distribution. 

r 

Subject to the above assumptions, the resulting "exact" equations 
are solved by an approximate method using a simple iteration procedure. 
An integral for the a:dal velocity is derived in terms of Known upstream 
conditions and the downstream pressure. The pressure distribution is 
found by integrating the radial momentum equation and applying the 
energy equation at a boundary. The mass flow equation is used to solve 
for the constant obtained from the axial velocity integral. The iteration 
is started by assuming an axial velocity distribution and continued until 
the solution is reached, i. e. , until the velocity distribution produced by 
an iteration is the same as that assumed for that particular iteration. 
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II. FUNDAMENTAL SOUATIONS 

A. Development of Basic Differential Equations 

The differential equations are developed assuming steady, 
compressible flov/ of a fluid through an axially symmetric cylindrical 
channel of arbitrary shape. An infinite number of blades is assumed 
so that the velocity components are considered independent of angle 
about the axis of the machine. The fluid is assumed perfect with no 
losses of any sort except as these losses can be approximated by a 
’’polytropic efficiency” and the use of the polytropic exponent "n”. 

All boundary layer and viscosity effects are also neglected. 

In order to generalise the equations for an arbitrary (but 
reasonable) channel shape, it is necessary to consider a special type 
of blading. These blades, infinite in number, have no other effect on 
the flow than to turn it through a prescribed angle distribution. Thus 
the effect of downstream blades on upstream conditions is neglected. 
Attention is directed to the flow patterns Immediately before and after 
the blade rows for either the stator or the rotor. 

In Figure i , a stream annulus of small radial extent is shown, 
and absolute velocity components are considered at stations 1, 2, and 
3. Generally, there will be a shift of the stream surface in passing 
through either stator or rotor blading, and three components of 
velocity will be present at each station. 

With the above assumptions, the continuity equation is written 
in differential form as; 



e, ^ 



("a 



( 3 ). 



-9- 




FIG. I— CHANNEL AND STREAM ANNULUS 

For a rov 7 o£ stator blades, the energy equation must hold along 
a stream surface, and for a row of rotor blades, the energy equation 
must hold along a relative stream surface. Hence, in terms of 
absolute velocities at stations 1, 2, and 3, we have: 






(4) 
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(5). 



Neglecting the radial force, the radial momentum equation .nust 
hold at all stations. In absolute velocities and omitting subscripii.s, 
this equation is: 

/ dV _ _ // iM: ( 2 ). 

TJr~ ~ T dr '^dz 



If the flow is considered to be isentropic along a relative stream 
surface, the following relations tnust hold along this surface: 
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( 6 ) 

( 7 ). 



If a polytropic flow is considered, ^ will be replaced by n in 
Equations 6 and 7. The above equations are fundamental in the solution 
of the dii'ect an 1 inverse problems. 



B. Transformation to Non-dimensional For'm 



For ease in handling the equations, and to avoid the use of large 
numoers in calculation, the derivation is completed in non-dimensional 
form. The notation used is given in Appendix A. Axial and radial co- 
ordinates are made non-dimensional by referring them to a standard 
dimension, usually the outer radius of the annulus at station i. Velo- 
cities are made dimensionless by referring them to a standard tip speed. 
The following notation is used; 

- outer radius of the annulus at station 1. 

- reference blade speed. • 

- dimensionless radial component of velocity. 

- dimensionless tangential component of velocity. 

- dimensionless axial component of velocity. 
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- dimensionless pressure. 

- dimensionless density. 

- dimensionless temperature. 

- gas equation. 
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- dimensionless mass flow. 



Wcj 




- dimensionless radial distance . 
•> dimensionless axial distance. 



With this notation, Equations (3), (4), (5), (2), (6), and (7) become 
respectively: 






( 8 ) 






J_^ 

Q 






dO 

3 ^ 




Qs _ 
Qa 




( 11 ) 



( 12 ) 



(13). 



C . Equa tions for th e Direct P ro blem for an Arbitrary Channel Shape 

The basic integral for the axial velocity will be derived for the 
rotor (stations 2 to 3) for. the isentropic case and simplified to the 
integral for the stator . The complete derivation is given in Appendix 
B for the isentropic direct problem and in Appendix C for the polytropic 
direct problem. 

Differentiating the left side of Equation (10) in respect to the 
radius ^ and the right side in respect to ^ gives: 
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(14). 



Using Equations (8) and (13), differentiation of the first term of the right 
side of the above equations results in: 
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r~! 
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(15). 



Substituting Equations (11) and (15) into (14), and using the differential 

cL 4& 



form of the nriass flow relation (Equation (8) ) for the ratio 

\ 

the following differential equation iB obtained: 
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(16). 



For a cascade with an infinite number *of blades, the ejcit flow 
is parallel to the tangent of the camber line at the trailing edge, and 
its direction is independent of the direction of flow approaching the 
cascade. For most practical applications, the direction of flow 
leaving the cascade is very nearly constant and independent of the 
inlet direction throughout the design operating range. ith this 
assumption for the rotor, it is therefore assumed that the rotor 
spouting angle jij is a known function of radius at station 3. Using 
the sign conv’ention shown in Figure 8, it is therefore possible to 
describe the relation between and along the radius at 

station 3 with a function O ( fj) such that: 
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Elirn.inating from the differential Equation (1 6) by use of 

Equation (17) and integrating, the following integral for the rotor axial 
velocity for the direct problem is obtained: 






rfj 






Ij: d.(4x.) 















(18). 



It is to be noted that no linearisation has been performed in 
arriving at the above integral. The expression for 4:^ is quite compli- 
cated, but its numerical evaluation by a simple iteration process to be 
described later is straightforv^ard. Once a set of streamlines have 
been assumed, all factors Involving the upstream conditions become 
known functions ol .In the solution of the direct problem, it 

is assumed that G is a given function of . Thus in the 

cL & 3 

integrand, P-, and ~r~i — unknowns which must be deterinined 

o oc. <3 

in the iteration process in order to evaluate the integral. 

An alternate form for the square bracket in the integrand can be 
obtained by using the energy and radial momentum equations at station 2: 

09). 



2^ -S ! -L A- (x^+d2\ - 224*^2)- JLt 

In the case that = 1 (isentropic expansion to station 2 from stagna- 
tion conditions), the bracketed term in the integrand may be further 
simplified to: 






( 20 ). 
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The constant of integration Cj is determined by integrating 
Equation (8). The integrated form of this relation equates the total 
flow rate upstream to that at station 3: 




( 21 ). 



Using Equation (13), the integral form of the radial momentum 
equation becomes: 




( 22 ). 



The constant of integration is obtained by applying the energy 
relation (Equation (10) ) at a boundary. 

Having assumed a set of streamsurfaces for a particular 

dL 6 

boundary configuration, the value of ^ can be approximated 

from the values of the slope and rate of change of slope of these stream - 
surfaces at station 3. At each point on a streamline at station 3, the 
follov.;ing relation must hold; 



v/herc 
the value of 
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(23), 



is the slope of the streamline in question. Thus 



dLGs 



becomes: 



cL0s ^ d4>. 
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In order to approximate the value of 
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(24). 



at points across the 



channel, it is therefore necessary to estimate the slope and curvature 



of the assumed streamlines at each point. 
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For the solution of the direct problem for the rotor v;ith an 
arbitrary channel shape, Equations (10), (13), (17), (18), (21), (22), 
and (24) are applicable. 

The derivation of the integral for the axial velocity in the case 
of the stator is carried out in a similar fashion qsing the form of the 
energy relation given as Equation (9), Using the sign convention of 
Figure 8, Equation (17) for the stator (stations 1 to 2) becomes: 

(^z) <Pa (25). 



The integral for obviously takes the same form as 
Equation (13), and is given here only for completeness; 



j-P^ ^cL4z 











■ r dZNr 


<p,f, [ 


T-, cLl 


, / 


■V 

dLQe 








(26). 



For the solution of the direct problem for the stator with an 
arbitrary channel shape. Equations (9), (12), (21), (22), (24), (25), 
and (26) (with appropriate subscripts) are applicable. 



^ • l^q^ations for the Inverse Problem for an Arbitrary Channel Shape 

In the Inverse problem for the rotor, the same assumptions and 
basic differential equations are used as in the direct problem with the 
exception that now the blade relation. Equation (17), is replaced by 
the following given function which describes the work o'ltput of the 
rotor as a function of the radius at station 3: 



( 27 ). 



- 16 - 



The function H gives thf; cimensionless rate of change 

of angular momentunn. per unit mass flov; at the particular radius, and 
represents the energy taken out of the stream. By multiplying this 
function by suitable constants it can readily be converted to B. T. U. 
per pound of mass flow. 

The complete derivation of the axial velocity integral is given 
in Appendix D for the isentropic case and in Appendix E for the poly- 
tropic process. For the isentropic inverse problem, the derivation 
remains the same up to and including Equation (16), which may be 
written in a slightly modified form as: 



r-t 






(28). 



Using Equation (27) to eliminate 



, the integral is 



found to be: 
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(29). 



Again the numerical evaluation of this integral by the iteration 
process to be described later is straightforward once a set of stream- 



ed 



appear in 



lines lias been assunoed. The unknown P_ and . 

3 

the integrand as in the direct problem. Because of the nature of the 
differential equation, hov/ever, the axial velocity < 2 also appears in 



the integrand. 
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£incc the mass flow and radial momentum integrals remain 
unchanged. Equations (10), (13), (21), (22), (24), (27), and (29) are 
applicable to the solution of the isentropic inverse problem. 

If the inverse problem could be conceivably applied to a stator 
it would consist of prescribing values of as a function of ^ , 

and solving for the resulting velocity field, distribution of energy and 
blade shape (as described by the blade trailing edge angle). In this 
case, the integral for the axial velocity simplifies to: 




(MJjl 



d&2 









(30) 



where is a prescribed function of the radius at station 2. 

Equations (9), (12), (21), (22), (24), and (30) (with appropriate 
subscripts) would be applicable. 



E. Changes in Equations for the Compressor Solution 

In the development of the equations in the preceding sections, 
the turbine sign convention given in Figure 8 has been adopted. If 
the equations arc to be used for compressor design. Equation (17) for 
the direct rotor problem becomes, in accordance vnth Figure 9: 

I 

K- &(g *s 

and the energy relation, Equation (10), is changed to: 
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Using these equations , the derivation is carried out in exactly the 
same manner as set forth in the preceding sections for the turbine. 

For the inverse problem for the compressor, the work out- 
put function (Equation (27) ) becomes: 

I 

Using Equations (32) and (33) in place of Equations (10) and (27), the 
relations for the inverse compressor problem are derived in the sanne 
manner as for the turbine. 
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in. PROPOSED ITERATION PROCEDURE 



A, Procedure for the Direct Problem 



The general scheme of iteration is exactly the same for the 
direct and inverse problems, subsonic or supersonic, isentropic or 
polytropic. In Appendix F, one complete iteration is given for an 
isentropic supersonic solution (Example Vb) and a polytropic subsonic 
solution (Example VIII b) for the direct rotor problem. The stator solu- 
tion is obtained by the same procedure using the appropriate axial 
velocity integral and energy equation. 

In order to understand the technique of iteration for the direct 
problem described below, the reader is urged to refer frequently to 
Examples Vb and VIII b of Appendix F, For an arbitrary channel 
shape, seven equations are used in the solution. For the rotor, these 
relations are Equations (10), (13), (17), (18), (21), (22), and (24). In 
the examples -worked out in Appendix F, the functions F and G 
have been assumed to be constant and independent of radius in order 
to simplify the axial velocity integral and the computation. This 
assumption is by no means necessary. For arbitrary given functions, 

F and O , the exponential integrals of Equation (18) 

should be evaluated at various stations across the channel by plotting 
the integrand and determining the value of the integral graphically. 

In the solution of the direct problem, it is assumed that the 
boundary configuration, blade speed, blade trailing edge angle, and 
upstream conditions rare prescribed. It is further assumed that the 
mass flow rate is fixed and is compatible v;ith the given upstream 



conditions. 
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Oncc the problem has beci set up, the first step is to assume 
a set o£ streamsurfaces as shown in Figures 3, 4, and 5. For most 
applications, these streamsurfaces should be chosen on the basis of 
equal annular areas in a plane perpendicular to the axis of the machine. 
The accuracy of the solution can be improved by Increasing the number 
of stations across the channel. However, it is believed that for most 
purposes, six or seven steps will give sufficient accuracy for design 
work. 

Having assurr'.cd a reasonable set of streamlines, the second 
step is to estimate the axial velocity distribution at the downstream 
station. Taking average or middle -channel values for the prescribed 
upstream conditions, the mass flow relation and the energy equation 
can he combined into one equation in terms of •i'j- Using average 
values in Equations (10) and (2i), the relation is obtained as follows: 






(35). 



Substituting Equation (34) into Equation (35) gives: 



m 






(36). 



Using Equation (J.7) to eli.r.inate A.? and Equation (23) to 
eliminate , Equation (36) becomes a relation involving only 

one unknown, (p^ . The solution must be found by trial and error. 

In general, there will be two values of <p^ which satisfy this 
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equation corresponding to the average subsonic and supersonic solutions. 



mass flow cannot be forced through the channel because of choking, 
and the boundary configuration or upstream conditions must be altered 
accordingly . 

Having found the average value of the axial velocity at station 3 , 

an estimate of the distribution of this velocity across the channel can 

be made considering the total temperature gradient, the blade shape, 

and the pressure distribution expected at the downstream station. 

Using the estimated distribution of , Equation (22) is 

iterated across the channel to determine the pressure gradient. The 

tangential velocity is found from Equation (17), the radial velocity from 

oL ^ 

Equation (23) , and the value of ^ from Equation (24) . In order 




0 — ^ 



d> 



FIG. 2 -PLOT OF MASS FLOW VS. AVERAGE AXIAL 
VELOCITY FOR CHANNEL FLOW 



If no value of ^ Datisfies Equation (36), the prescribed 



to evaluate ^ 




oLO^ 




at various points across the 



, and 
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chann@l, it is necessary to estimate as accurately as possible the slope 
and curvature of each streamline at the station in question. Knowing 
the complete channel configuration, slopes and curvatures at points 
across the channel can be easily estimated using a simple linear distri- 
bution scheme, or a large scale drawing can be constructed and the 

values of slope and curvature taken graphically from this. The value 
cL / V 

of - j -f, ' in Equation (24) can be taken as the average value between 

CL fs 

stations along the axis. However, it is to be noted that even a rela- 
tively large error in the estimation of the above quantities will have little 
effect upon the pressure gradient obtained because of the nature of the 
integrated form of Equation (22). The actual pressure distribution 
across the channel is found by applying the energy equation at a boundary 



, the axial velocity 



to determine the constant € 2 * 

With these values of P, and 

3 aLCs 

integral can now be iterated across the channel for values of <{>2 in 

terms of the constant C If the curvature of the channel is not large, 
aL03 

is very small compared to the other factors in 



the value of 






the integrand ( 2 G, for example), and hence a relatively large error in 
cLO^ 

estimating . has a small effect on the resulting value of 9 -. 

The constant C . is obtained by iterating the mass flow relation 
(Equation (21) ) across the channel, and a new ajdal velocity distribution 
is obtained. 

The above described process is the procedure for one complete 
iteration of the given equations. Normally in the subsonic case, the 
distribution of <>52 thus obtained is used in the next iteration, and tlic 
process is repeated until the axial velocity distribution produced by an 
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iteratioa is the same as that assumed for that particular iteration. 

This is the solution. Values of and are then found by 

use of Equation (13). 

B. i^^rocedure for the Inverse Problem 

Generally, the same procedure of iteration, as described in 
section A above for the direct problem, is used in the solution of the 
inverse problem. Equations (lO), (13), (21), (22), (24), (27), and (29) 
are applicable to the isentroplc rotor. Relations used in the poly- 
tropic case are derived in Appendix E. Examples IX, X, and XI of 
Appendix F are typical solutions to the inverse problen . One complete 
iteration is given in Example IX b to illustrate the method used. 

It is noted that the integrand of Equation (29) contains a term 
involving . Since in any practical solution the axial velocity 

at any point across the channel will not be aero, this term does not 
cause any real difficulty in the iteration procedure. However, if the 
original estimate of the axial velocity profile at the downstream station 
differs greatly from the actual solution, it has been found that succes- 
sive iterations will "oscillate" about the real solution, but will 
eventually converge as in the direct problem. 

It is obvious that the work output function H {^ 3 ) » given as 

Equation (27), must be prescribed reasonably. For any given channel 
shape, blade speed, upstream conditions and mass flow, there is a 
limited range of blade angles that can be used before the smooth flow 
over the blades is destroyed by separation effects. The limitation thus 
imposed by the maximum turning angle through the blade row also 
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limits not only the amount of work that can be extracted, but also the 
distribution of this v./ork across the channel. 

C . Convergence of Subsonic and Supersonic Solutions 

Figure 2 indicates that, for any given setup, there is normally 
a supersonic as well as a subsonic solution for the direct and inverse 
problems. This is illustrated in Appendix F by Examples IV a and V a 
for the stator and Examples IV b and Vb for the rotor. In these 
examples, for a given boundary configuration, blade spouting angle, 
and set of upstream conditions, both the subsonic and supersonic 
solutions have been obtained. 

Experience in v/orklng various examples has shown that if the 
flow is in the low subsonic range, point A of Figure 2, for example, 
the iteration process converges very rapidly, and only two or three 
iterations are necessary to obtain accuracy of four significant figures 
in the solution. However as the flow becomes highly subsonic as at 
point B, the convergence becomes less rapid, and more Iterations are 
necessary to obtain the solution. It is interesting to note that, in the 
subsonic range, the process will converge from any reasonable 
assumed distribution of even possible to obtain a solution in 

' which the flow Is supersonic in one part of the channel and subsonic in 
the remainder (Examples la and III a). 

It has been found, hov/ever, that if the flow is completely 
supersonic (point E of Figure 2, for example), the iteration process 
diverges rapidly away from the solution. Thus if an axial velocity 
distribution is assumed whose average value is represented by point D, 
successive iterations will result in decreasing values of and 
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eventually the subsonic solution at point A will be reached. If an 
axial velocity distribution corresponding to point F is assumed, 
successive iterations will rapidly increase 4*3 toward infinity. It 
is the realization of this phenomenon which makes it possible to 
obtain the isentrcpic supersonic solution in a given problem, fiaving 
made two successive iterations at each of the points D and F, and 
having noted the relative change in the constant in these iterations, 
a simple interpolation will give the values of <J>^ corresponding to the 
solution at point E once the solution has been "straddled". 

It is important to understand the basic reason why the proposed 
iteration process is successful. For either the direct or inverse 
problems, six unknowns must be solved for at the downstream station. 
These are the three components of velocity ( ^3 \ ^3 ^md the 
pressure P^, the density * and the temperature ^ . Using 

the isentropic (or pol^-tropic) relation along streamsurfaces, two of 
these (Qj and 2 ^ ) can be eliminated. The remaining four unknowns 
are made functions of the axial velocity and the boundary and 
streamline configuration. The radial velocity is eliminated in terms 
of <j>2 by Equation (23), the tangential velocity by Equation (17) 

or (27), and the pressure by the integral. Equation (22). Thus, 
effectively, only one unlcnown remains, the axial velocity (j)^ , and for 
any assumed distribution of all other unknowns are determined. 

An iteration process must be used to determine 93 because of the 
complexity of the equations Involved. 
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D. Correction for the Disolaccment of Streamlines 



Having assumed a set of streamsurfaces for any particular 
example, and having found the corresponding solution by the proposed 
iteration procedure, it is possible to make a correction to determine 
more accurately the actual position of the streamsurfaces and thus 
refine the solution. The mass flow functions m f and ™ (fs) 
are plotted against radius for conditions at stations Z and 3 as shown 
in Figure 6. Using this diagram and fixing the radial coordinates at 
station 3, it is possible to estimate new corresponding upstream 
coordinates. V/ith these corrected upstream radial coordinates, and 
the upstream conditions corresponding to them, the problem is re- 
worked and a new solution found. It is to be emphasized that this is 
a small correction and that the nev.^ solution will not vary more than 
one or two percent from the solution based on the original stream - 
surfaces. An Illustration of this method is given in Example V b of 
Appendix F. 
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IV. ADDITIONAL CORRECTIONS AMD SI^-.IPLIFICA TIOMS 

A. Corrections to the Basic Solution 

If it is desired to solve for the flow pattern and distribution of 
energy in the field at any axial station in the blade row itself, this can 
be done approximately by multiplying the axial velocity by a function 
which accounts for the blade thickness and thus reduces the effective 
area of the channel. Assuming a blade thickness function t(^) and 
a blade spacing s , the increased axial velocity at any point witliin the 
blade row would be: 



position across the channel. The assumptions used in the direct 
problem, including that of axial symmetry, remain unchanged, except 
that the blade spouting angle now becomes the mean camber line at 
the point in question. The fundamental equations and iteration pro- 
cedure remain unchanged except for the use of the above indicated 
expression for the axial velocity. 

One of the basic assumptions used in the derivation of the 
fundamental equations given in Part II, is that the effects of the vorti- 
city of the flow dov;nstream on upstream conditions are negligible. A 
correction to the upstream conditions due to vorticity effects down- 
stream could be made using the approximate theory suggested by 
Marble (Reference i5) for the three-dimensional case. If the hub 
ratio is large, the simpler two-dimensional approximation given by 
Rannie (Reference 1) can be used. Once the upstream conditions have 




(37) 



where s and t (C) are assumed to be known functions at any radial 
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been changed using this higher order correction, the corresponding 
new solution can then be found by the same procedure as previously 
outlined in Part III. 



B. Simplification of Flow Patterns 

An investigation has been made to determine whether or not 
there are any special cases in which the equations, used in this analysis 
of three-dimensional, compressible flow, reduce to a simple expression 
for the axial velocity. Consider the integrated radial momentum 
equation: 




Since the value of appears in the integrated mass flow 
relation. Equation (21), to the power -y , it is clear that, generally, 
no simplification can be made for a channel of arbitrary shape where 



the values of 






d&3 






and 






dd3 

d^- 



cannot be expressed in 



terms of constants and the radius ^ 

However, one simplification does occur in the example of a 

rotor in a cylindrical channel. In this case, if the velocity adter the 

r 

rotor is purely axial. Equation (22) reduces to in 

addition, the total temperature and pressure are constant across the 

channel at station 2, and if (vortex flow). Equation 

7Z 

(i6) for a straight channel reduces to ~ Cj. The constant C -2 is 
evaluated by applying the energy equation at a boundary, and the mass 
flow integral becomes: 
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been changed using this higher order correction, the corresponding 
new solution can then be found by the same procedure as previously 
outlined in Part III. 

Simplification of Flow Patterns 

An Investigation has been made to determine whether or not 
there are any special cases in which the equations, used in this analysis 
of three-dimensional, compressible flow, reduce to a simple expression 
for the axial velocity. Consider the integrated radial momentum 
equation; 




Since the value of appears in the integrated mass flov.' 

relation, Zquation (2i), to the pov/er -y , it is clear that, generally, 

no simplification can be made for a channel of arbitrary shape where 

the values of 0^ cannot be expressed in 

terms of constants end the radius ^ 

However, one simplification does occur in the example of a 

rotor in a cylindrical channel. In this case, if the velocity after the 

r 

rotor is purely axial, Equation (22) reduces to = O-^-i . If, in 

addition, the total temperature and pressure are constant across the 
channel at station 2, and if ^ (vortex flow). Equation 

72 

(l6) for a straight channel reduces to cfjj = Cj. The constant C -2 is 
evaluated by applying the energy equation at a boundary, and the mass 
flov^ integral becomes: 



-29- 



m 



I 




(38) 



or, since Tzt is constant 




(39). 



Even in this simplified case. Equation (39) must be solved by 
trial and error for the constant C,. This particular example gives a 
constant mass flow per unit area and the work is taken out uniformly 
across the channel. 

With the exception of the above case, it can be generally stated 
that there are no special conditions or flow patterns in three-dimen- 
sional, compressible flow through turbomachinery using this method 
which result in a simple closed-form solution for the downstream axial 
velocity. 

C, Proposed Optimuirt Turbine Design Condition 

In the design of a turbine, the condition v;hich normally controls 
the amount of heat that can be released in the combustion cliamber is the 
limiting temperature that the first stage rotor blade root can safely 
withstand in continuous operation. Since the blade stresses decrease 
from root to tip, the allowable temperature is higher at the tip than at 
the root. This fact allows the possibility of increasing the amount of 
heat released in the combustion chamber, and hence the energy in the 
jet exhaust, for a fixed amount of work taken out of a turbine. 



In formulating the optimum turbine design condition, reference 
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will be made to a one -stage turbine v/lth an arbitrary, but fixed, 
boundary configuration, such as is shov,/n in Figure 3. It is assumed 
that the total pressure in the combustion chamber, the blade speed, 
mass flov/ and total amount of work to be extracted by the turbine are 
fixed by conditions in the compressor. It is further assumed that the 
distribution of the maximum allowable total temperature across the 
channel at station 2, based on limiting blade temperatures, is known. 
Since the total enthalpy is conserved through the stator row, the limiting 
total temperature profile at station 1 is therefore known. 

In order that the maximum energy be left in the exhaust gases 
at station 3, the condition is imposed that the velocity shall be axial 
after the turbine ( = 0) and that the total enthalpy shall be constant 

across the channel at this station ( ” ®)- 

cLfj 

To simplify the formulation of the proposed optimum turbine 
design condition, it is assumed that the known limiting total temperature 
profile at station 1 is linear, and that the co-nbustion chamber can be 
constructed so as to give this desired total temperature variation across 
the cliannel. 

The design problem can then be formulated as follows: Given a 
one -stage turbine with fixed boundary configuration, combustion chamber 
total pressure, blade speed, rpass flow, and linear total temperature 
profile at station I, what are the blade spouting angles of stator and 
rotor so that a given amount of work can be extracted by the rotor and 
yet leave the total temperature constant and the velocity axial at station 
3? 

With the above assumptions and conditions. Equation (27) can be 



written: 
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^(0 = 



(40) 



and Equation (lO) in terms of total temperatures becomes: 




If ^3f is to be constant at station 3, then the work must be taken 



out across the channel in such a way that: 




(42). 



Since the linear change in total temperature at station 2 is 
assumed known, the distribution of work across the channel is deter- 
mined by Equation (42) and the total amount of work to be extracted. 
The solution to this simplified inverse problem is found in a manner 
similar to that Illustrated by Example X of Appendix F, and the blade 
trailing edge angles of stator and rotor arc completely determined. 

In view of the discussion given in Part III, Section B, on 
limiting blade angles, it will obviously not alv;ays be possible to obtain 
the above proposed optimum condition, particularly if the upstream 
total temperature gradient is large. Furthermore, the above consi- 
derations could not reasonably apply to a multi-stage turbine since 
boundary layer effects and mixing would soon destroy the original 
total temperature profile. However, the designer should not overlook 
the possibility of increasing the gross thrust of a turbine by the use of 
a total entlialpy gradient based on blade limiting temperatures. (See 
Example Xlf Appendix F. ) 
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V. resuijTS and conclusions 

An approximate method has been developed for the solution of 
three-dimensional, compressible flow in axial turbomachines. The 
methods and equations used are extremely flexible and can be used to 
solve the direct and inverse problems, subsonic or supersonic, isen- 
tropic or polytropic for either a turbine or a compressor. The upstream 
conditions, blade slmpe (as defined by the blade trailing edge angle), 
blade speed, and axial boundary configuration can be prescribed in the 
direct problem, while the upstream conditions, blade speed, blade 
loading, and axial boundary configuration can be prescribed in the 
inverse problem. In each case, six unknov/ns are solved for at the 
downstream station. These are the three components of velocity, the 
pressure, the density, and the temperature. In the analysis, maximum 
utilization has been made of the given channel configuration. 

The equations are derived for a blade row with infinitely many 
blades, on the basis of isentropic (or polytropic), non-viscous, com- 
pressible, axially- symmetric, steady flow. Radial or nearly radial 
blades are assumed so that the radial force can be neglected. Separa- 
tion of flow and boundary layer effects are not considered in the analysis. 
For an arbitrary boundary shape, a special type of blading is assumed 
whose only effect on the flow is to turn it through a specified blade 
angle distribution across the channel. Using the assumption that the 
blade shape is specified by the blade trailing edge angle only, and 
assuming that this spouting angle is very nearly independent of condi- 
tions ahead of the blade row in the operating range, the analysis of the 
flow has been made at the station immediately dov^nstream of the blade 



- 33 - 



row in, terms of specified conditions immediately upstream of the blade 
row. 

This method greatly simplifies the analysis of flov.’ in turbo- 
machines, since, with the above assumptions, it is not necessary to 
solve for the complete three-dimensional flov/ field through the blade 
row, but only to consider stations immediately upstream and down- 
stream of the blading. 

The iteration process developed is simple and straightforward, 
and the same generrJ. procedure is used in the solution of all suggested 
types of problems. The calculations can be done with a slide rule, and 
the accuracy of the solution can be improved by increasing the number 
of radial steps across the channel. Corrections can be made for the 
displacement of the assumed streamlines and for the effect of down- 
stream vorticity on upstream conditions to further improve the solution. 
It is believed that the results thus obtained by this approximate method 
are sufficiently accurate for design purposes. 

Once the boundary configuration and other necessary conditions 
have been assuiricd for any particular problem, experience in working 
various examples by the proposed iteration process has shown that a 
basic solution for the direct problem accurate to four significant 
figures can be obtained for a stator in approximately four hours and 
obtained for a rotor in approximately six hours. The calculations for 
inverse problem arc somewhat simpler and less time is normally 
required for this solution. 

cLO 

In order to evaluate the term cL ^ ' appears in both 

the integral for the axial velocity and the integral for the pressure in 
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the case of an arbitrary channel shape, it is necessary to estimate the 
slope and curvature of the streamlines at the station in question. It is 
to be noted that, for a reasonable channel shape, a relatively large error 
in this estimation will have little effect on the solution obtained. The 

de 



term 



dL^ 



appears in the general axial velocity integral for the 



rotor (Equation (IS) ) and for the stator (Equation (26) ) as a small 
correction term, and a large error in the assumed values of slope and 
curvature will have a small effect on the resulting axial velocity profile. 
This is especially true for the rotor where the term 2G is normally 
dominant. 

Considering the pressure integral (Equation (22) ), it is to be 

dd 



noted that here, again, the contribution of the term 






is small 



for a reasonable boundary shape. For the stator, the term 



is 



the predominant one, and the terms 



o dO j ^ <LO 

9 -tt and <P -tt are 

usually small corrections. Because of the small magnitude of the factor 



multiplying the integrand of Equation (22) in the iteration, a relatively 

dd 



large error in 



dLC 



will have but a slight effect on the resulting 



pressure gradient. After the rotor, where the velocity is predominantly 

axial, the term might be of the same order of magnitude as 

H ft d^Q 

9 ^ " and (p . As a result, it is possible to have a 

df cL^ 

negative pressure gradient radially after a rotor if the slope and curva- 
ture of the channel is large. 

It is impossible to make any general statement as to the relative 
contribution of the integral (Equation (18) ) and the constant of integration 
Cj to the resulting axial velocity profile. The relative size of the 
integral and the constant may be of any magnitude depending on the 
channel shape, total temperature gradient, blade spouting angles, and 
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apstream conditions of any particular example. 

Several interesting observations are possible from the results 
of the examples given in j^vppendix F. The channel shapes assun.ed for 
these iterated problems are given as Figures 3, 4, and 5. In each 
example, the assumptions, channel configuration, upstream conditions 
radial coordinates, and other data are specified, and the”resulting 
solution at the downstream station is tabulated. One complete iteration 
is given in Examples Vb, VIII b, and IX b to show the method used in 
each case. 

Examples I and II use the channel configuration shown in 
Figure 3. In both of these examples of the direct problem, the mass 
flow, blade speed, blade spouting angles, and total heat content in the 
flow are the same. Example I, however, has no radial total tempera- 
ture gradient at station 1 , while Example II has a large radial enthalpy 
gradient. The effect of a total temperature gradient across the channel 
in changing the axial velocity profile Is clearly evident from the results 
obtained. Normally for a stator, it is to be expected that the positive 
pressure gradient would make the axial velocity higher at the root than 
at the tip. This is the case in Example I a for no temperature gradient. 
However, Example II a shows that the effect of the positive total 
gradient imposed is to make the axial velocity higher at the tip than at 
the root, and reduce the Mach number slightly at all stations across the 
channel. Also, the resulting total pressure at station 3 has been in- 
creased by the use of the radial enthalpy gradient. 

Examples IV and V use the channel configuration of Figure 5, 
and have the same mass flow, blade speed, blade spouting angles, and 
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conditions at station 1, Example IVa is the subsonic stator solution 
v;hile Example Va is the supersonic solution for the same stator. 
Examples IVb and Vb are the subsonic and supersonic rotor solutions, 
respectively, for conditions after the subsonic stator of Example IVa. 

The mass flow functions for the subsonic and supersonic rotor examples 
are plotted in Figure 6, From this graph, the position of the displaced 
streamlines has been estimated and plotted on Figure 7, This plot of 
the streamline shift in each case, clearly indicates that the shift is in 
one direction when the ’'governing'' Mach number is subsonic, and in the 
opposite direction when the "governing" Mach number is supersonic. 

The "governing" Mach number has been defined by Monroe (Reference 7, 
page 51) as the Mach number based on the meridional velocity instead 
of the total velocity. The results shown in Figure 7 clearly agree with 
Monroe's observation that the deflection of the " ixaeridlonal Mach 
surfaces" is in different directions for subsonic and supersonic 
"governing" Mach numbers, and that the deflection of the streamsurfaces 
will be less if the governing velocity is transonic than if it is entirely 
subsonic or entirely supersonic. 

Examples VI, VII, and VIII use the channel configuration of 
Figure 5, the same conditions at station I, and the same blade speed, 
blade spouting angles, and mass flow. Example VI is the isentropic 
case. Example VII is a polytropic solution with n = 1. 37 (^ = 94. 5 per- 
cent), and Example VIII is a polytropic solution with n = 1. 33 
( y = 86. 8 percent). From these examples, the effect of a decrease 
in "polytropic efficiency" on the three-dimensional velocity field and 
distribution of energy in the flow can clearly be seen. When the 
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" efficiency” has dropped to 94. 5 percent, the average axial velocity at 
station 3 has increased 2. 7 percent, the average total pressure has 
decreased 2. 5 percent, and the work output has increased 4. 7 percent 
when compared with the isentropic solution. As the "polytropic 
efficiency” drops to 86. 8 percent, the results of Example VIII show that 
the average axial velocity at station 3 has increased 8. 2 percent, the 
average total pressure has decreased 8. 0 percent, and the work output 
has increased 17.0 percent when compared with the isentropic solution. 

The results of these examples clearly show the reason for the 
"loading up” of the first stages of a multi-stage turbine as the efficiency 
drops. As the efficiency decreases, more work is taken out in the 
initial stage, and less total pressure is left for the remaining turbine 
stages. Thus, the last stage of a multi-stage turbine may well act as 
a compressor if the efficiency of the system has dropped sufficiently. 

In view of this phenomenon, it is obvious that the design engineer should 
design his machine for the expected efficiency in the operating range, 
and not on the basis of isentropic flow. 

Example XI has been constructed to illustrate the increase in 
gross thrust obtainable by the use of a radial total temperature gradient. 
Assuming a one -stage turbine configuration similar to that shown in 
Figure 3, and keeping the mass flow, combustion chamber total pressure, 
blade speed, and work extracted constant in each case, the gross thrust 
of the jet exhaust has been calculated with and without a total tempera- 
ture gradient. In order to give the example more physical meaning, 
blades were used which were tapered in cross-sectional area from hub 
to tip. Having calculated the centrifugal stress at various radial points 
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along the blades, and using the high temperature data for ''Vitalliam” 
(n. S, 21) (NR -10) given in References 16 and 17, it v/as estimated 
that the rotor blade root could safely stand a temperature of 1450°F, 
■wMle the blade tip temperature could be 1850*^F. 

Using this data, the gross thrust of the turbine was calculated 
for: (a) A constant radial total temperature of 1450°F, and (b) a linear 
total temperature gradient of 400°F with the same root total tempera- 
ture of 1450°F. The results of these examples show that an increase 
in gross thrust of 917, 8 pounds, or 7. 7 percent, was obtained by the 
use of the radial total temperature gradient indicated. 
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APPENDIX A 



SYMBOLS AND DEFINITIONS 



c 

Cl 



a(fj) 



Stator spouting angle measured from the plane perpendicular 
to the axis. 

Rotor spouting angle measured relative to the rotating wheel. 
Axial distance between stations in feet. 

Constant obtained from axial velocity integral. 

Constant obtained from radial momentum integral. 

Specific heat at constant pressure. 

Used to denote incremental changes of a quantity in the radial 
direction. 

Cotangent of stator spouting angle ^9^ as a function of radius. 

Cotangent of rotor relative spouting angle a.s a. function 
of radius. ^3 + 4s 



Y 

g 

m 

M 

Mg 

n 

7 

w 

p 

p 



Ratio of specific heats. 

Acceleration of gravity (32. 17 ft/sec^). 

Dimensionless worK output of rotor as a function of radius. 



Dimensionless mass flow. 



m - 



zn^p^Ro 



$ 



- Mach number based on total velocity. M - 

- "Governing" Mach number as defined on p. 87 of Reference 7. 

- Polytropic exponent. 

- Polytropic efficiency. ^^2/ ~ ri^i 

- Angular velocity of rotor in radians per second. 

- Pressure in pounds per square foot. 

y> 

- Dimensionless pressure. P = 



‘.3 



P 

Q 

R 



T 

T 

u 

9 

V 

w 

4» 

I 

V/ 

z 

< >1 
< >2 

< >3 

< >i 

< )o 

Bf 

d^-f 



7 4 

- Density (lb. sec. ‘‘lit. ). 

- Dimensionless density. Q ~ 



7°t 



- Gas constant (1715 



ft^ 



). 



- 2o~ 

sec F 



- Reference radius in feet. 

- Radius from the axis in feet. 

~ Dimensionless radius. ^ 

- Absolute temperature in degrees Rankine. 



R 



T 



- Dimensionless temperature. ^ 

- Radial component of absolute velocity in feet per second. 



u. 



- Dimensionless radial component of absolute velocity. 9 - ' 

- Tangential component of absolute velocity in feet per second. 

- Dimensionless tangential component of absolute velocity. A - ^ ^ 

- Axial component of absolute velocity in feet per second. 

- Dimensionless axial component of absolute velocity. <j> - (J ^ 



- Dimensionless total absolute velocity. § = * 9^ 

- Mass flow in pounds per second. 

- Axial distance along axis of rotor from reference point in feet. 

- Dimensionless axial distance, t 



- Subscript to denote station immediately upstream of stator. 

- Subscript to denote station immediately downstream of stator 
and immediately upstream of rotor. 

- Subscript to denote station immediately after rotor. 

- Subscript to denote inner boundary of channel. 

- Subscript to denote outer boundary of channel. 

- Slope of streamline. 

- Rate of change of slope of streamline. 



ti 
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APPENDIX B 



Note: In Appendices B, C, D, E and F, Equations of lower number 
than 43 refer to relations previously given in the text. The 
numbering of equations is consistent throughout the thesis. 
Development of Equations for Direct Problem (Isentropic). 

The equations are derived for a rotor with infinitely many blades 
on the basis of isentropic, non-viscous, compressible, axially-symmet- 
ric, steady flow. Radial or nearly radial blades are assum.ed so that 
the radial force can be neglected. The channel shape, blade shape, 
mass flow, and upstream conditions are assumed specified in any par- 
ticular example. 

The derivation follows the general procedure given in Refer- 
ence 1, Section 2:2. 2 for the incompressible case. The energy equa- 
tion is v/ritten for any two points (for example, and b^. Fig. 5) on 
a streamline surface between stations 2 and 3. Using absolute veloc- 
ities, the equation in dimensionless form is: 




where 




and 




The energy equation is now differentiated in respect to the rad- 



ius at stations 2 and 3: 



45 









Using the gas equation, Z'j = , the first term of the right 

side of Equation (44) becomes, with the iscntropic relation: 



- (aX - 

4^2 \Pal ■ [raJ 



(13) 



■y clTs _ _y_ ) _ _2l ^ / Ay p ^ ) 

y-, " /'-/ cL^y Qsi y-‘ dL^3 \Qa j 



(45) 



Using the differential form of the continuity equation for a 
stream cylinder of small radial extent, 

( 8 ) 



and differentiating Equation (45) as a product, the follov/ing is obtained: 



-2L ^ 

y-, d^, 







-21 

r-i 



cL 





Using the radial momentum equation at station 3, 



J_ _ ?\3 _ n ddj _ j d &3 



(15). 



( 11 ) 



and combining Equations (8), (11), (15) and (44), the following differ- 
ential equation results: 




1 
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J-^ + ^ ^ jL eL dilii2^3) 

Qs ^fs fo ^ \ ^a I ^&4e ^ ^ J ^ ^3 



f y drat _ dii^z) 

V-' ■ 









i !^ -(b ^ 

^3 



^- /\^ , , d(^s)\s) 



= ( ^-3^ <Psf 3 [ _2l £{^2^ JL £ 3 , 3. ^ (pH 

^ Pj ~ ^-> Qs P^r d4^\ Q^ J 



<^fs 



(46) 



For the direct problem it is assumed that the rotor spouting 
angle is known as a function of radius at station 3. It is there- 

fore possible to describe the relation between Aj and (ps along 
the radius with a function G(^) such that: 



^3 - G^s) i>3~ ^3 



(17) 



Eliminating from Equation (46) with this relation the fol- 
lowing differential equation is obtained: 



iff, 



_ ^ PsJ^ Ps$3 






r'-t u. 

JL _il t(P^\ ^ i—(^ P^ ) 






dOa 

^3 



G%3 



H>s .£ + 



.(PjY ^3 r r atar d(4a)ia) (n Pah 






dP>s . ^ [_Gf_ / dOf-a P) 

^^3 ^('-^Vafs 



-m 



Psf fs 



y dTat 



a( 4a 
P^a 



X afj^Q'J 



I , 3^3 i o /O 

5-77 i-^6 



(47) 



m 




t 
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If P,, 4, j and the knovm quantities at station 2 are 

considered to be functions of (to be determined in the iteration 
process) the integral for <f>j can be immediately written down as: 
43 






dfs 

^ ^ 



r 1 


/ 


/ r 


l<- 





y cLZat 



'T-/ 






L-\ 






dLOa 



d4. 



'■) 






(18). 



An alternate form for the square bracket in the above inte- 
grand can be obtained by using the energy and radial momentum equa- 
tions at station 2: 








(19). 



In the case that p 2 ^ = 1 (isentropic expansion to station 2 from 
stagnation conditions), the bracketed term in the integrand may be 
further simplified to: 



[ 






( 20 ). 



The integral form of the mass flov/ equation is: 

and the integral form of the radial momentum equation is: 



( 21 ), 
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Equations (10), (13). (17), (18), (21), (22), and (24) are used 
in the solution of the direct problem for the rotor (isentropic). 

In the case of a stator. Equation (17) becomes: 






(25). 



G is replaced by F in Equation (18) and the terms 



and 



cL (4z ) 

2G are omitted. Equation (18) becomes Equation (26) for the isentropic 



stator. Equations (13), (21), (22), and (24) (with appropriate sub* 
scripts) remain unchanged, and the energy relation becomes Equation 
(9). 



APPENDIX C 



Development of Equations for Direct Problem (Polytropic). 

The equations for a polytropic process through a rotor are 
derived with the same assumptions as used in Appendix B with the 
exception that the isentropic relation (13) is now replaced by: 



J. / 




where n is related to the so>called polytropic efficiency rj for an 
expansion (n < y) by: (see Reference 12, p, 449) 



^77 = (49). 

72-/ / y-/ 

The energy equation (10), the continuity equation (8), the rad- 
ial momentum equation (11), and the blade relation (17) remain the 
same in differential form. 

The first term of the right side of Equation (44) .now becomes; 



_2l SL- (-3.) - _21 _L 

(PsCifj 

y cL f Pj j _ / cL Ps , / / d. Rj _ y Ps dQs 

y-t Qjd-^3 , 



But, 



/ / cLRs _ y Ps dQs 

y-t y-r d^^ 



ri/ 



—P 



r d 




So, as before, 



i 

/ 

] 

( 



I 
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The right side of the differential eqixation is obtained in the 
same form as Equation (46) with a correction term as follows: 



> j 




(A) - 



r-i 
y p r 



8 



7^ 

Yn 



I Q^J 



Pa d 

Q. df^ 



71- t " 

P 



(51). 



Using the differential form of the continuity equation (8), Equation 
(51) finally becomes: 



7--/ -5 




r-/ 's d^^ )\ Ps / 4>^ (r-o (^- 0 [<pzJdLf/ 



(52). 



The differential equation corresponding to Equation (46) for 
the isentropic case now is: 



X ± —{"X^ + d>^\ +■ J_il ilA _ 



^4 



_ /' 

y-/ ^ 



3T'] 

\ Qa ) 









^(r-p(72-/) \ a fy 



m 



yt-i 



a P3 



(53). 



The relation corresponding to Equation (47) becomes: 



^ +<p 



B(H-G^) 






_ / Pa '^ ^3 r y dTzt _ 9(4 z X^ 
\PeJ y-/ 






_JLp^d— (At) 

y-/ ^ V / 









(7r-^)(7i-/) 



(4) 



Hzi 
dS ^ 



( 54 ). 



The integral equation corresponding to Equation (18) is: 



i: 






II 









<0 




51 






f’ d(, 

o 



■ ft 



! 4a 

\|77^ \Pzl <Pzf& 



■y' cLZet _ ol-(^sX ^ 



/^) + J. -^j: 2IL(S£)^^ df-f-C 



(55). 



If = 1 at station 2, the square bracket under the integral 



will be: 



r cLTst f / _ dJ^eXt) 

r-> d.4,(‘ J cL4^ 



(56). 



The integrated form of the mass flow equation is 
'%■ 

and the integral form of the radial momentum equation is: 

A 



Ps ‘ 



Tl-/ ( Qa \ / ( n cL03 j. ci^J \ ^ j ^ 

^ \¥lj (■^ 

— 4t 



_, _ 2 _ 
n-i 



(57) 



(58). 



Equations (10), (17). (24), (48), (55), (57) and (58) are used 
in the solution of the polytropic direct problem for the rotor. 

In the case of the stator. Equation (25) is used in place of 
Equation (17), and the terms -rr (^z Xg') and 2G are omitted in 

GL 

Equation (55). Equations (10), (24), (48), (57) and (58) remain un- 
changed (with appropriate subscripts). 
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APPENDIX D 



Development o£ Eciuations for the Inverse Problem (Isentropic). 

The following equations are derived with the same assump- 
tions as given in Appendix B with the exception that now the blade 
relation. Equation (17), is replaced by tiic following given function 
which describes the work output of the rotor as a function of the 
radius ^3 : 



For the isentropic case, the derivation remains unchanged 
up to an including Equation (46), which may be written as: 




(27). 






hs , ! d 





(28). 



Using Equation (27) to eliminate , the differential eqixa 



tion corresponding to (47) is developed as follows: 



^3 



hs . / cL X 3 4s^3 cL(43^4) 

j:. • ^ "■ jt s t , 











( 59 ) 




1 
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Assuming that the right side of the above Equation (60) is 
a known function oi f 3 (to be determined in the iteration process), 
the integral equation for tl\e asial velocity can be written as: 

. / (PsY 4s ^ r dz,t r p^d P^^\ A V . 

i- 

(29). 



Since the mass flow integral and the radial momentum inte- 
gral remain the same. Equations (10), (13), (21), (22), (24), (27) and 
(29) are used in the solution of the isentropic inverse problem for 



the rotor. 



I 



I 



\ 

I 

I 

I 



1 



I 
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APPENDIX E 

Developr>-^ent of Equations for the Inverse Problem (Polytropic). 

The following derivation is made assuming the same conditions 
for a rotor as given in Appendb: C. The work output of the rotor as a 
function of radius is again given by Equation (27). 

The differential equation (53) still holds, which is: 

2 






r~, 's ^(r-OCn-j) \ Qjdfp J 



(53). 



The resulting integral equation is similar to Equation (29) with 



the addition of a correction jterm: 



( PsT is J oLTzt _ _2L n^'cL f_^\ , i \ 















( 61 ). 



The mass flow and radial momentum equations remain the 
same as in Appendix C, and Equations (10), (24), (27), (48), (57), 
(58) and (61) are used to solve the polytropic inverse problem for 



the rotor. 



Ill 




I 
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APPENDIX F 
Iterated Examples 

The following data are used for Examples I-X inclusive: 



= 734 rad/sec (7000 RFIv'".) = 16,300 



lb 



ft‘ 



R = 1.2500 ft, 
o 



Y 

R 



= 1 .40 



= 1715 



fr 

sec^°F 



(j R^ = 917.5 ft/sec 



cj^PJ = .8418 J£ 10^ ~ 



C = .2396 



BTU 

lb°F 



c = .250 ft 



sec 



Example la . 

Solution for the stator of Fig. 3 (cylindrical channel), with no 
radial total temperature distribution and coordinates for equal areas 
in the plane perpendicular to the axis . This is a solution for the direct 
problem ) given) for the isentropic, subsonic case. 

For the stator with F constant, = I, and Q = 0, Equations (18), 
(21), (22), (13), (25) and (43) of Appendix B become: 



= 



/ 






% O^F^) 









c, 






fp 



(62) 






TTL = 









(63) 



Pz 



?<-/ / 



^4 ^ 
^2 



r> 



f± 






% 






r 



(64) 



f • m- It) 



Tz l r-/ 



(65) 
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/^2 = (F = const.) 

^ (V.ach number) 

Conditions at station 1; (given) 



( 66 ) 



(67) 



( 68 ). 





P,J= 1.00 


^"^It ~ 


.23828 


Tn,= 4.1968 


4>J = .4857 




r-, = .9723 

JL 


^1 = 


.23354 


r, = 4.1631 


A/ = 0 


dZt 


0 F = 2.75 


m 


= ,020417 


W = 114.7 


T,. = 2060 
sec It 


Colution at station 2; 












^2 


^2 


^2 


*^2 


i. 




.8000 


.8544 


.9055 


.9539 


1 .0000 




.7926 


.7479 


.7105 


.6785 


.6508 




2.1796 


2.0567 


i .9539 


1.8659 


1 .7897 


^2 


.4928 


.5365 


.5730 


.6041 


.6307 




3.4283 


3.5128 


3.5794 


3.6340 


3.6739 


^2 


.14372 


.15274 


.16007 


,16623 


.17144 


N-' 


1 .059 


.986 


.929 


.880 


.841 




+.0001 


- .00C2‘ 


-.0001 


-.0004 


-.0002 


-1 = 


.6508 


TD - p 

‘ 2t " -^It 




c = c 

2t It 


®2 = 0 



C, = . 



8169 




I I 
i I 






I 

i 

I 
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Note: The column headed t,X ^ above represents the difference 

between ^ calculated by the iteration process and 
obtained by applying the energy equation along the stream- 
lines indicated after the solution has been reached. The same 
applies in the solution of the following examples . This is a 
rough check on the accuracy of the solution. Obviously, if an 
exact solution had been found, A T would be zero. 

Example Ib . 

Solution for the rotor of Fig. 3 (cylindrical channel) with no 
radial total temperature gradient and radial coordinates for equal areas. 
This is an isentropic, subsonic solution for the direct problem (G {^3 ) 
given) for a rotor. 



For the rotor with G constant and = 1 , the applicable 
relations are Equations (63), (64), (65) (with appropriate subscripts) 
and the following three equations: 




~ ^ (G = const.) 



(70) 

(71) . 



Conditions at station 2 are the same as after the stator , Example la , 



cL 'Czt 



= 0 



G = ! .411 



m = .020417 



I 

I 
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Solution at station 3: 





^3 


^3 


^3 


^3 


®3 




. 8000 


. 8544 


.9055 


.9539 


1. 0000 


<P3 


. 6624 


.6779 


.6955 


. 7146 


.7346 


}\3 


. 1346 


. 1021 


.0759 


. 0544 


.0365 


^3 


. 5862 


. 5865 


. 5865 


. 5865 


. 5865 




3.6025 


3.6034 


3.6034 


3.6035 


3.6035 


<^3 


. 16269 


. 16278 


. 16276 


. 16276 


. 16277 


^3t 


.6244 


.6254 


.6264 


.6273 


.6284 




3.6679 


3.6699 


3.6717 


3.6734 


3.6751 


^3t 


. 17019 


. 17041 


. 17061 


. 17077 


. 17099 




0 


-.0007 


-.0017 


-.0035 


-.0058 



Cj = . 5710 C 2 = . 8585 

Eaiample 11a. 



Solution for the stator of Fig. 3 (cylindrical channel) with the 



same 


conditions and assumptions 


as given in 


Problem la except with a 


radial total temperature gradient. 


1 




~ 




Equations (62) - ( 68 ) of Example la are applicable. 




Conditions at station 1: (given) 










^1 


^1 




^1 


®1 


c 


.8000 


. 8544 


, 9065 


.9539 


1.000 




3.3319 


3.8252 


4. 2415 


4.6361 


5.0118 




l.OO 


1.00 


1.00 


1.00 


1.00 


^it 


.29569 


.26142 


. 23577 


.21570 


. 19953 


J’l 


.9723 


.9723 


.9723 


.9723 


. 9723 


Xt 


3. 3548 


3.7946 


4.2076 


4. 5990 


4.9717 


Qj 


. 28981 


. 25622 


.23103 


.21141 


. 19556 


<t>, 


.4355 


.4628 


.4871 


. 5096 


. 529C 




0 


0 


0 


0 


0 




F = 2. 75 


m = .020417 


= + 8 . 1495 
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Solution at station 2: 





“2 


^2 


*^2 




®2 


^3 


. sooo 


. 0544 


. 9055 


.9539 


1. 0000 


4>z 


. 6778 


.6805 


.6829 


.6848 


.6864 




1.8639 


1.8714 


1. 8780 


1. 3832 


1.8876 


^2 


. 5295 


. 5713 


. 6061 


.6355 


.6611 




2. 8200 


3.2599 


3.6762 


4.0729 


4. 4527 


®2 


. 18777 


. 17525 


. 16488 ' 


. 15602 


. 14845 


M 


.999 


.930 


.883 


.842 


.805 




0 


-.0012 


-.0052 


-.0104 


-.0172 






C j = . 5565 


= .8339 








^2t = ^It 




»®2t ~ ^It * 


02 = 0 





Example lib 

Solution for the rotor of Fig. 3 (cylindrical channel) with the 



•*ame conditions and assumptions given in Problem Ib except with a 
radial total temperature gradient. Equations (63), (64), (65), (69), (70) 
and (71) of Example I (with appropriate subscripts) are applicable. Condi- 



tions at station 2 are the same as after the stator, Ei^ample lla, above. 





G = 1.411, 


cLTzt 


- 8. 1495, 


m a .020417 




Solution at station 3: 












®3 °3 






S 


®3 


6 


. 8000 . 8544 




.9055 


. 9539 


1.0000 




. 5377 .6174 




.6920 


.7628 


.8303 




-.0413 .0168 




.0709 


. 1224 


. 1716 


^3 


.6009 .6009 




.6009 


.6009 


.6012 


^3 


2.9237 3. 3072 




3.6670 


4.0081 


4.3334 




.20551 .18168 




. 16386 


. 14990 


. 13871 


^3t 


.6313 , .6381 




.6422 


.6447 


.6461 


^3t 


.21289 .18964 




. 17184 


. 15759 


. 14603 


Tat 


2.9653 3.3643 






4.0895 


4.4235 




0 +.0035 




+.0012 


-.0039 


-.0126 




Cj = .4635 




C 2 = .8646 
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After correcting the upstream conditions at station 2 for the 
displacement of the streamlines from the position originally assumed, 
the following corrected solution at station 3 is obtained: (See Problem 
Vb for method. ) 



<3 


^3 




*^3 


^3 




®3 


8000 


. 8544 


.9055 


.9539 


1. 


0000 




. 5311 


.6106 


. 6853 


. 7560 




8234 


^3 


-.0506 


0072 


0615 


1128 


. 


1618 


P 3 


. 6029 


6029 


. 6029 


. 6029 




6032 




2.9268 


3. 2771 


3.6296 


3.9868 


4. 


3374 




. 20602 


. 18411 


. 16623 


. 15127 


. 


13917 


^3t 


.6327 


. 6378 


.6417 


.6453 




6476 


'^3t 


2. 9675 


3. 3303 


3.6951 


4 . 0648 


4. 


4263 


^3t 


. 21325 


19166 


. 17381 


15879 


• 


14641 




0 


-. 0001 
Cj = .4578 


0022 
^2 = 


- 0055 
. 8654 


•*' • 


0113 



Example Ilia - 

Solution for the stator of Fig. 4 (divergent channel with straight 
walls) with a radial total temperature gradient. This is an isentropic, 
subsonic solution for the direct problem for a stator (F ) given). 

For the stator with F constant and = 1, Equations (63) , (65) 
and (66) are applicable in addition to the following four relations: 




t 6i 



r 




Conditions at station i: (given) 









^1 




®1 




,80 


, 85 


.90 


.95 


1.00 


T^it 


3.9931 


4.0950 


4. 1968 


4. 2987 


4. 4006 


®lt 


. 25043 


. 24420 • 


.23828 


. 23263 


.22724 




l.OO 


1.00 


1.00 


1.00 


1.00 




.9696 


.9696 


.9696 


.9696 


.9696 


r, 


3.9581 


4.0591 


4. 1600 


4.2610 


4. 3620 


Qi 


. 24497 


.23888 


.23309 


.22756 


. 22229 


<^1 


. 4944 


. 5014 


. 5069 


. 5112 


, 5142 


0. 


-. 0247 


.0000 


.0253 


.0511 


.0771 




0 


0 


0 


0 


0 




F = 2.75 


s+2.0375 m = 


.024639 V/ = 


13s!4 



c = . 25 ft 






i 



I 
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Solution at station 2: 





^2 


*^2 


*^2 


^2 


^^2 


. 793 


850 


• 910 


• 970 


1. 30 




. 7410 


. 7046 


. 6735 


. 6469 


. 6236 


'^2 


2. 0377 


1 9376 


1 8521 


1. 7790 


1. 7149 


9 ^ 


-.0370 


0 


. 0337 


. 0647 


. 0935 


^2 


. 5247 


. 5716 


. 6122 


. 6472 


. 6777 




3. 3212 


3.4904 


3. 6479 


3. 7966 


3. 9376 




. 15798 


. 16378 


. 16732 


. 17051 


, 17210 


M 


1. 004 


.930 


. 870 


. 817 


.779 




0 


0027 


0061 


0103 


-. 0138 




Cj = . 6017 


= .8317 


P = 
2t 


^It ' ^2t "" 


oi., = 



Example Illb . 

Solution for the rotor of Fig- 4 (divergent channel with straight 
walls) with a radial total temperature gradient. This is an isentropic, 
subsonic solution for the direct problem for the rotor (G ( given). 
Equations (63), (65), (7C) , (73) , (75) (with appropriate sub- 



scripts) and the following two relations are applicable: 




Conditions at station 2 are the same as after the stator, Example 
Ilia above. 

dZzt: 



G = 1.411 m = .024639 



= 1. 6979 c = . 25 ft. 



i 



I 



I 

I 
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Solution at station 3: 



^3 


^3 
. 780 


^3 
. 850 


*^3 

.920 


^3 
. 990 


®3 

1. 060 




. 4567 


.4987 


. 5423 


. 5873 


. 6347 


03 


-. 0228 


0 


. 0271 


. 0587 


. 0952 




-. 1356 


-. 1463 


- 1548 


-. 1613 


-. 1644 


^3 


.6500 


. 6500 


. 6505 


. 6508 


. 650S 




. 18409 


. 17954 


17525 


. 17119 


. 16720 


rs 


3. 5308 


3. 6209 


3, 7117 


3. 8027 


3. 8923 


^3t 


. 6712 


.6749 


. 6780 


. 6804 


.6828 


Qst 


. 18836 


. 18442 


. 18051 


. 17672 


. 17304 




3. 5634 


3. 6600 


3.7560 


3. 8514 


3.9459 




+. 0001 


+. 0005 


-. 0013 


- 0048 


-. 0092 



Cj = . 3871 = . 8842 

Example IVa 



Solution for the stator of Fig. 5 (divergent curved channel wails) 
with, a radial total temperature gradient. This is an iaentropic, sub - 
sonic solution for the direct problem for a stator (F given). The 

same equations as used in Example Ula are applicable. 



Conditions at station 1: (given) 
»1 ‘■l ■=! 



f/ 


. 80 


. 84 


.88 




3. 7894 


3. 9524 


4. 1152 


Pit 


1. 00 


O 

o 


1. 00 


Oit 


. 26389 


. 25301 


. 24299 


Pi 


. 9400 


.9400 


.9400 


r, 


3.7231 


3. 8832 


4. 0433 


^1 


. 25248 


. 24207 


. 23248 


4>i 


. 6804 


. 6960 


. 7090 


®i 


-. 0340 


0 


. 0354 




0 


0 


0 


F 


= 2. 75 


If 


4. 0746 
c = . 



.92 

4. 2783 
1. 00 
. 23374 
.9400 
4. 2034 
. 22363 
.7205 
. 0720 
0 



.96 

4.4413 
1.00 
. 22516 
.9400 
4. 3636 
. 21542 
. 7293 
. 1094 
0 



1. 00 
4. 6043 
1. 00 
. 21719 
.9400 
4. 5237 
. 20780 
. 7365 
. 1473 
0 



m = . 029242 V.- 



164.2 



sec 



c = . 250 ft. 
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Solution at station 2: 





^2 


»^2 


^2 


^2 


®2 


^2 




,785 


.842 


.899 


.956 


1 .013 


1 .070 




.7051 




.6531 


.6344 


.6192 


.6072 




1 .9390 


1.8598 1 


.7960 


1 .7446 


1.7028 


1 .6698 


«2 


-.0705 


.0135 


.0914 


.1649 


.2353 


.3036 


^^2 


.5417 


. 5860 


.6222 


.6518 


.6761 


.6962 




3.1805 


3.3928 3 


.5937 


3.7860 


3.9713 


4.1519 




.17032 


.17272 


.17313 


.17217 


.17022 


.16767 




.977 


.909 


.854 


.809 


.776 


.747 




0 


0 

Cj = .5694 
- I- 

2t It ’ 


.0013 
^3. = 


-.0039 -.0069 

= .8393 

’ ^2t ■ ^It 


-.0118 



Example IVb 

Solution for the rotor of Fig. 5 (divergent curved channel walls) 
with a radial total temperature gradient. This is an isentropic, 
subsonic solution of the direct problem for the rotor (G ^ ) given) . 
The same equations as used in Example Illb are applicable. 

Conditions at station 2 are the same as after the stator , 
•Example IVa . 



♦ 



dn 



at 






G = 1.730 



= 2.8593 



m 



029242 
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Solution at station 3: 
^3 *^3 


^3 


^3 


®3 


^3 


6 


.760 


. 848 


.936 


1.024 


1. 112 


1. 200 


*3 


. 1958 


. 2641 


. 3260 


. 3852 


.4458 


. 5117 


^3 


-.4213 


-. 3911 


-.3720 


-.3576 


-. 3408 


-. 3148 




0294 


. 0106 


. 0750 


. 1617 


. 2719 


.4094 


P 3 


. 6948 


. 6981 


. 7010 


. 7027 


.7010 


. 6927 




3.4149 


3.5669 


3.7180 


3.8682 


4. 0126 


4. 1461 


^3 


. 20346 


. 19573 


. 18852 


. 18167 


. 17466 


16707 




3. 4660 


3. 5998 


3.7535 


3.9064 


4. 0568 


4. 2018 


^3t 


. 21117 


. 20027 


. 19303 


. 18619 


. 17950 


. 17272 


^3t 


. 7320 


.7209 


.7246 


. 7274 


. 7283 


. 7257 




. 212 


. 210 


. 219 


. 236 


- 266 


. 302 




0 


+. 0010 


-. 0003 


0049 


-. 0113 


-. 0196 



Cj = . 1594 C 2 = .9012 

Example Va - 

Supersonic isentropic solution for the stator of Example IVa 
(divergent curved channel walls and radial total temperature gradient). 
The conditions at station 1, the data, and tlie applicable equations are 
the same as for Example IVa. 



Solution at station 2: 



^2 


^2 


^2 




®2 


^2 


02 .9303 


.8906 


. 8601 


.8371 


.8201 ^ 


3087 


^2 2. 5583 


2.4491 


2. 3653 


2. 3020 


2. 2553 2. 


2239 


-.0930 


. 0178 


. 1204 


. 2176 


.3116 


4043 


.3171 


. 3715 


.4169 


.4547 


.4856 


5109 


2. 7294 


2.9784 


3 2051 


3.4157 


3.6131 3. 


7990 


Q 2 .11617 


. 12471 


. 13007 


. 13310 


. 1 3440 


13430 


M 1.367 


1. 274 


1. 188 


1. 123 


1.080 1. 


038 


0 

0 


+. 0039 


+. 0034 


-. 0012 


-.0082 


0178 


Cj = . 7512 


C 2 = J203 


^2t ■ ^It' *^2t ■ ^ 


If 





I 

i 
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Example Vb . 

Supersonic isentropic solution for the rotor of Example IVb 
(divergent curved channel walls and radial total temperature gradient) . 

The conditions at station 2 are assumed the same as after the subsonic 
stator of Example IVa. The data and equations used are the same as that 
used for Example IVb. 

One complete approximation for this problem is given below to 
show the method of iteration. The procedure of correcting for the dis- 
placerrent of the assumed streamlines is then given to show the method 
for this correction. 

The iteration is started by assuming an axial velocity distribu- 
tion across the channel. The tangential velocity is then calculated using 
Equation (70). The value of ^ ^ is found from Equation (75) and 
the pressure distribution from Equation (73) and the application of Equa- 
tion (77) at the inner channel boundary. The axial velocity distribution 
is then determined using Equation (76) and the constant is determined by 
the mass flow Equation (63) . The iteration is continued until the values 
of <j>j assumed agree with those resulting from the iteration process. It 
is important to note that in the supersonic case* the iteration diverges 
rapidly away from the solution. Thus if an axial velocity distribution is 
assumed that is greater than the solution, successive iterations will in- 
crease 4*2 imtil the pressure has reached aero. If an axial velocity 
distribution is assumed which is less than the solution, successive 
iterations will result in decreasing values of and eventually the subsonic 
solution will be reached. The supersonic solution is found by a simple 
interpolation once the solution has been "straddled". 
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The reverse is true for the subsonic solution for a given setup. 
The iteration process will converge toward the solution from any 
reasonable assumed distribution of 

Data used in the solution: 

“y cLVet 

,2 



= 



088 



= 10. 0075 G= 1.730 2G= 3.460 






1>G‘ 



= . 7496 



^3 

(If) 

Factor 

/ 

^•7496 

i fc-Af,) 

— i 



^3 


’=3 


‘^3 


^3 


®3 




.760 


.848 


.936 


1. 024 


1. 112 


1. 200 


-. 150 


.040 


.230 


.420 


.610 


. 800 


-.250 


. 100 


.450 


. 800 


1. 150 


1. 500 



r ^ 



02117 


.02125 .02133 


.02138 


.02144 


4310 


1.5415 1.6411 


1.7305 


1.8098 


i 1. 


> 

1316 1.0509 .9824 


9235 .8722 


01299 


.01357 .01489 


.01694 


.01719 


7702 


-.8526 -.9333 - 


1.0018 


-1. 0842 


006498 


.006235 .005993 


.005768 


.005560 



2t 



(**Factor*’ above = 



Try Cj = 1 . 39 20 



5 



4 



1 : 






4 









I 
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H 


^3 


*^3 


d3 


*^3 




<Ps 


1.7098 


1. 6727 


1 . 6943 


1.7698 


1.9004 


2.0912 


G4>, 


2.9580 


2.8938 


2.9311 


3. 0618 


3. 2877 


3. 6178 




2.1980 


2.0458 


1.9951 


2. 0378 


2. 1757 


2.4178 


&3 


-. 2565 


.0669 


.3897 


.7433 


1.1592 


1.6730 


a 


4.8312 


4. 1853 


3.9804 


4. 1526 


4.7337 


5. 8458 




6.3568 


4.9354 


4.2525 


4. 0553 


4.2569 


4.8715 



e, 

303 

-'g 



5. 6461 


4. 5939 


4.1539 


4. 1561 


4. 5642 


-. 0948 


. 2283 


. 5665 


.9512 


1.4161 


3.6750 


3.6681 


4.0181 


4.7261 


5.8386 


+. 3484 


-.8374 


-2.2763 


-4.4955 


-8, 2680 



d4>3 


5.0235 


4.9820 


5. 2060 


5.6770 


6.4060 


7.4200 




-.7535 


. 1993 


1. 1974 


2. 3843 


3.9077 


5.9360 




-,AZ1A 


.1673 


.7624 


1.4158 


2. 1855 


3. 1368 


d03 


-1. 1809 


.3666 


1.9598 


3.8001 


6.0932 


9.0728 


-4>^ 


f2.0191 


-.6132 


-3.3205 


-6.7254 


-11. 5795 


-18.9730 



- 4 >^ 



f. 7030 -1.9668 -5,0229 -9.1524 -15.2762 



7 "-/ 



^ K 04352 + .01116 - .01885 - ,05475 - .10553 
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The constant C2 is found by applying Equation (77) at the inner 
channel boundary; 

= 3.7894-. 14286 [4. 831 2+2. 9234f. 0658+2(1. 5221 H. 6705)] = 1.7600 
P3 = . 5417(. 5534)T?7 = .0683 = (.0683)^^= .4645 



r-y 




Average 



^3 

.4645 


^3 

.5080 


*^3 

■ .5192 


^3 

. 5003 


®3 

. .4456 


'3 

.3400 


.0683 


.0934 


. 1008 


.0386 


.0591 


.0229 


. 1261 


.1594 


. 1620 


. 1359 


. 0874 


.0329 


.2278 


.2694 


. 2725 


.2404 


. 1754 


.0873 


• 


2486 


2709 


2564 


2079 


1313 


. 5355 


.4920 


.4309 


.4996 


. 5543 


.6601 



.5137 .4864 .4902 .5269 .6072 



The axial velocity integral is iterated in steps across the channel as 
follows : 



a 




Average 
. 04877 


. 09754 


.09754 


.16448 


.13388 


.23142 


.31687 


. 17090 


.40232 


. 50862 


.21260 


.61492 


. 74602 


. 26220 


.87712 





The mass flow equation is iterated to solve for Cj as follows: 

. 029242=. 003229(. 04877+1. 1799Cj)+. 003676(. 16448+1. 091 2Cj) 

00381S(. 31687 +1. 0166Cj)+. 003522(. 50862+. 9529Cj) 

+. 002257(.74602+. 8978Cj) 

— .023790 I 

7 mm ~ 



I 

!_ I 

) 




I 
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Solution at station 3: 





®3 


^3 


*^3 




®3 


^3 




1.7098 


1. 6727 


1.6943 


1.7698 


1.9004 


2. 0912 


As 


2. 1980 


2. 0458 


1.9951 


2. 0378 


2. 1757 


2.4178 


O 3 


-.2565 


. 0669 


. 3897 


.7433 


1. 1592 


1. 6730 


^3 


. 0683 


. 0934 


. 1008 


. 0886 


.0591 


. 0229 




1.7601 


2. 0079 


2. 1365 


2. 1406 


1.9793 


1. 5653 


^3 


. 03880 


. 04653 


.04718 


.04139 


. 02986 


. 01464 


Ut 




3. 0093 


3. 1204 


3.2056 


3.2572 


3. 2649 


^3t 


. 13256 


. 12794 


. 12162 


.11358 


. 10373 


. 09198 


^3t 


.3815 


. 3849 


.3795 


.3641 


.3379 


. 3000 


M 


1.780 


1. 570 


1. 540 


1.620 


1.880 


2.420 




1. 100 


.993 


1. 005 


1 . no 


1. 340 


1.810 


KOOOl 


f. 0028 


-. 0170 


-.0548 


-. 1065 


-. 1603 



Correction of supersonic solution for displacement of streamlines: 



II 








^2 ^2 


^2 




®2 


'2 


^2 . 785 . 842 


.899 


.956 


1. 013 


1. 070 


0 .005498 


.011201 


.017071 


. 023095 


. 029242 


%■ 


^3^3^ f 3 








aj bj 

. 760 . 848 


^3 

.936 


*^3 

1.024 


®3 

1. 112 


^3 

1. 200 


?n.(^^) 0 .005258 


. 011609 


.018374 


. 024639 


.029242 



The functions m(^) and ) are plotted against radius on 



Fig. 6. From this graph, a corrected ^ is obtained and, hence, 
new conditions at station 2. 
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Corrected values at station 2: 





^2 


^2 


*=2 




®2 


'2 


4a 


.785 


.839 


.902 


.968 


1. 026 


1. 070 


4>z 


.7051 


.6778 


.6322 


. 6305 


.6157 


. 6072 




1.9390 


1 . 8640 


1.7926 


1.7338 


1. 6933 


1.6698 


9^ 


-.0705 


. 0090 


.0955 


. 1804 


.2514 


. 3036 




3.7894 


3. 9438 


4. 1240 


4.3127 


4.4785 


4. 6043 


^2 


. 5417 


.5837 


.6241 


.6580 


.6816 


.6962 


tz 


3. 1805 


3. 3815 


3. 6043 


3.8266 


4.0135 


4. 1519 




. 17032 


. 17259 


. 17315 


.17197 


. 16978 


. 16767 


Corrected solution at station 3: 










^3 


•^3 


^3 


*^3 


®3 




^3 


1.7398 


1.7007 


1.7151 


1.7686 


1.8722 


2.0493 


^3 


2. 2499 


2.0942 


2.0311 


2.0357 


2. 1269 


2. 3453 


03 


-. 2610 


. 1020 


. 2573 


.4422 


.9361 


1.6394 


P3 


.0605 


.0842 


.0953 


.0959 


.0733 


.0295 




1.7003 


1.9450 


2. 1067 


2.2072 


2. 1223 


1.6832 


O3 


. 03558 


. 04330 


. 04523 


.04344 


.03452 


.01754 




2. 8660 


2.9896 


3. 1189 


3.2376 


3. 3064 


3. 2898 


®3t 


. 13104 


. 12665 


. 12076 


. 11325 


. 10449 


.09375 


^3t 


.3643 


.3790 


.3763 


.3666 


.3458 


.3071 




1-.0006 


+.0033 


-. 0069 


-.0368 


-.0881 


-. 1032 



Cj = 1.4165 C 2 = .4487 

Note: Examples VI, Vll and VUl use the channel configuration of Fig. 
5 and the same conditions at station 1. Example VI is the isen* 
tropic solution (y = 1.40), Example VII is a polytropic solution 
for n = 1. 37, and Example VIII is for n = 1. 33. From these 
examples it is possible to see clearly the effect of a decrease 
in polytropic efficiency on the solution for the same mass flow, 
upstream conditions and channel configuration. 






! 
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Example Via. 

Solution for the stator of Fig. 5 (divergent curved channel 
walls) with a radial total temperature gradient. This is an isentropic 
(v s 1. 40) subsonic solution for the direct problem for a stator 
(FC/g) given). The radial coordinates have been chosen for equal 
areas in the plane perpendicular to the axis. The same equations 
as used in Example Ilia are applicable. 



Conditions at station 1: (given) 





^1 


^1 


*^1 


‘'l 


®1 


'1 


4 , 


.8000 


.8438 


.8854 


.9252 


.9633 


1.0000 




3.7894 


3.9679 


4. 1374 


4.2995 


4. 4548 


4. 6043 


^It 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


Qlt 


. 26389 


.25202 


.24170 


. 23258 


. 22448 


.21719 


^1 


.9500 


.9500 


.9500 


.9500 


.9500 


.9500 


Oi 


. 25439 


.24295 


. 23300 


. 22421 


.21640 


.20937 


r, 


3.7343 


3.9102 


4.0772 


4. 2369 


4. 3900 


4.5373 


<P, 


.6202 


.6355 


.6481 


.6583 


. 6657 


. 6716 


9i 


-.0310 


.0030 


.0367 


.0701 


. 1026 


.1343 


F = 


2.75 m 


s . 026865 


dr,t . 


4. 0746 




0 


Solution at station 2: (7 = 


^ 1.40) 










^2 


^2 


*=2 


^2 


•=2 


h 


4 


.7850 


.8496 


.9098 


.9661 


1.0194 


1.0700 




.5772 


. 5513 


. 5317 


.5166 


. 5046 


.4948 




1. 5873 


1. 5161 


1.4622 


1.4206 


1. 3876 


1. 3607 


9z 


-.0577 


.0193 


.0865 


. 1453 


. 1986 


. 2474 


P 2 


.6712 


.7095 


.7378 


.7592 


.7763 


. 7898 


O 2 


. 19848 


. 19723 


. 19451 


. 19103 


. 18732 


. 18349 




3.3814 


3 . 5974 


3.7930 


3.9738 


4. 1437 


4. 3041 


^2 


.775 


.754 


.675 


. 668 


.630 


. 609 




0 


-.0013 


-. 0025 


-.0037 


-. 0061 


-. 0080 




It 


^2t ■ ^It 


' ^2t “ ^It ^ 


j = . 4661 


= .8923 



8923 
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Example VIb 

Solution for the rotor of Fig. 5 (divergent curved channel walls) 
with a radial total temperature gradient. This is an iscntropic (7=1.40) 
subsonic solution for the direct problem for the rotor (G(^) given). 

The radial coordinates are for equal areas. The same equations as 
used in Example Illb are applicable. 

Conditions at station 2 are the came as after the stator. Example 



Via. 

G=1.96 m = . 026865 4?^ = 2.8593 

Solution at station 3: (-y = 1.40) 





^3 




^3 


<^3 


®3 


^3 


^3 


.7600 


.8661 


.9605 


1 . 0464 


1.1258 


1. 2000 


4>s 


. 1652 


. 2370 


.2982 


.3543 


.4089 


.4658 


93 


•*. 0248 


.0187 


. 0844 


.1660 


. 2616 


. 3726 


^3 


-.4362 


-.4016 


-.3760 


-.3520 


-. 3244 


-.2870 




.7550 


.7595 


.7626 


.7635 


.7606 


.7532 


^3 


3.4970 


3. 6683 


3.8290 


3.9802 


4. 1197 


4. 2460 


O 3 


. 21589 


. 20707 


. 19916 


. 19181 


. 18460 


. 17738 




3. 5281 


3.6993 


3.8605 


4.0126 


4. 1549 


4. 2867 


^3t 


.7788 


.7824 


.7847 


.7853 


.7837 


.7788 


^3t 


. 22073 


. 21150 


. 20326 


.19572 


. 18859 


. 18167 


42} 


0001 


-. 0002 


-.0024 


-.0072 


-.0135 


-. 0219 



Cj = .1329 C2 = .9228 

Example Vila 

Solution for the same channel configuration, mass flow, and 
conditions at station 1 as for Example Via, except for a polytropic 
process (n = 1. 37 , y = 94. 5°/o) through the stator. 

~ Equations (48), (57), (58), (66), (74), (75) and the following 



integral are applicable for the polytropic stator (F = const): 



I 



1 





74 



t 






\T}^,\r-, cLf^i' )ycL^^ 



(78). 



, (r-?L) Tjt: <LPs ^ . /y.>c2 

(r-ocyi-o <P^ df 2 ^ oC^ 



F = 2.75 ms .026865 

Solution at station 2: (n = 1. 37) 



a 



cLT/j& 



= 4.0746 





^2 

. 7850 


^2 

.8496 


^2 

.9098 


^2 

.9661 


®2 

1.0194 


'2 

1.0700 




. 5897 


. 5652 


.5469 


. 5330 


. 5222 


. 5135 




1.6216 


1. 5543 


1 . 5040 


1.4657 


1.4360 


1.4121 




-.0590 


.0203 


.0890 


. 1499 


.2055 


.2567 




.6450 


.6837 


.7127 


.7347 


.7523 


.7662 


r^ 


3. 3635 


3. 5778 


3.7726 


3.9526 


4. 1218 


4.2814 


^2 


. 19176 


. 19110 


. 18889 


. 18585 


. 18251 


. 17895 


^2t 


.9789 


.9821 


.9845 


.9863 


.9874 


.9882 


^2t 


. 25834 


. 24751 


.23793 


.22938 


. 22164 


.21461 


aTz 


+.0001 


-.0017 


-.0022 


-.0038 


-.0066 


-.0090 




.4762 


C 2 * . 8883 


II 

44 







Example VHb 

Solution for the same channel configuration* blade shape* 
and mass flow as for Example VIb* except for a polytropic process 
through the rotor (n = 1. 37, ^ = 94. 5®/o). 

Conditions at station 2 arc the same as after the stator* 
Example Vila. 

Equations (48), (57), (58), (70). (75), (77) and the follow- 
ing integral are applicable for the polytropic rotor (G = const, 






I: 





I 

I 



75 










(r-n) 






- (A 

-/;l (^A 



72 V 



dPs^ 






G‘ 

/•>«-<j2 






i 



c 

g« 



G s 1.96 


m a . 026865 




2.8593 






Solution at station 3: 


(n = 1. 37) 












^3 




^3 


^3 


®3 


^3 




.7600 


.8661 


.9605 


1 . 0464 


1. 1258 


1.2000 


4>s 


. 1664 


.2496 


.3120 


.3637 


.4097 


.4514 


^3 


-.4339 


-. 3769 


3490 


-.3335 


-.3228 


-.3153 


Bs 


-.0258 


.0197 


.0883 


, 1704 


. 2621 


.3611 




.7387 


.7430 


.7457 


.7463 


.7439 


.7379 


Ts 


3.4890 


3.6590 


3.8190 


3.9694 


4. 1094 


4. 2382 


^3 


. 21172 


.20306 


. 19524 


. 18799 


. 18101 


. 17410 




3.5199 


3.6839 


3.8422 


3.9947 


4. 1404 


4. 2807 


^3t 


.7620 


.7608 


.7617 


.7632 


.7636 


.7640 


Qst 


.21646 


. 20653 


. 19823 


. 19102 


. 18443 


. 17849 


aTs 


0 


0044 


-.0092 


-.0137 


-.0177 


-.0194 


Cj = ..1338 




= .9215 









It is interesting to note that* compared to the isentropic case 
(Examples VIa,b), the polytropic flow through the identical stator 
and rotor results in the following percentage changes at station 3 
(maximum); 



‘>3- 


^5.3% 


T -0. 3**/o 


^ 3* 


+7. 7°/o 


-3.0% 


^3^ 


-2. 3°/o 


Workout: +6.6°/o 



1 




I 



m 
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Example Villa 

Solution for the same channel configuration, mass flow, and 



conditions at station 1 as for Example Via, except for a polytropic 
process (n = 1.33, TP = 86.8°/o) through the stator. The equations 
used in Example Vila are applicable. 

F = 2.75 m=. 026865 

Solution at station 2: (n ~ 1. 33) 



^2 

O, 



'^2t 



^2 


^2 




.7850 


.8496 


.9098 


.6225 


.5994 


. 5828 


1.7119 


1.6483 


1.6027 


-.0622 


.0216 


.0948 


. 5378 


.6280 


.6584 


3.3149 


3. 5286 


3.7225 


. 17731 


. 17796 


. 17685 


.9337 


.9470 


.9532 


. 24770 


. 23863 


.23035 


-.0001 


-. 0002 


-.0019 



Cj = . 5027 C 2 = 



'^2 


®2 




.9661 


1.0194 


1.0700 


. 5704 


. 5613 


. 5545 


1. 5686 


1. 5436 


1 . 5249 


. 1605 


. 2209 


.2772 


.6817 


.7000 


.7146 


3.9022 


4.0695 


4. 2279 


. 17470 


. 17199 


. 16902 


.9571 


.9608 


.9631 


. 22262 


. 21564 


.20918 


-. 0044 


-.0071 


-.0107 



8765 



Example Vlllb 

Solution for the same channel configuration, blade shape, and 
mass flow as for Example VIb, except for a polytropic process 
through the rotor (n = 1. 33, y = 86. 8°/o). 

Conditions at station 2 are the same as after the stator, Ex> 
ample Villa. The equations used in Example Vllb are applicable. 

One complete approximation for tlie polytropic, subsonic, direct prob- 
lem for the rotor is given below to demonstrate the iteration procedure. 
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Data used in the solution: 
m = .026865 G = 1 .96 



r dTzt 

r., 



= 10.0075 



2G = 3.9200 






1+G^ 



= .7935 






IH 






m 



Factor 



^3 




^=3 


<^3 


®3 


^3 


.7600 


.8661 


.9505 


1 .0464 


1.1258 


1.2000 


-.1500 


.0791 


.2829 


.4684 


.6398 


.8000 


-.2500 


.1720 


.5474 


.8891 


1 .2049 


1.5000 



r-, d( I 



fZ 

(r-yt) 

(r-/)(n 

/ 



L(3r) 



■793f 



q(i^£d 



.02083 


.01874 


.01716 


.01593 


.01495 


1 .6273 


1.7561 


1 .8551 


1.9330 


1.9949 


9.1313 


9.0490 


8.9897 


8.9174 


8.8593 


9.6600 


9.8311 


10,0646 


\ 

10.2042 


10.2714 


2.0529 


2.1450 


2.2328 


2.3169 


2.3975 


33 1.1208 1.0325 .9646 .9102 .8653 


.01532 


.01530 


.01516 


.01495 


.01470 


.006801 


.006099 


.005064 


.004509 


.004072 



(For definition of "Factor" above, see Problem Vb) . 



^ >a 

1 



I 
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Try Cj = 


.1568 






^3 


‘=3 


4>J 


.1949 


.2717 .3312 


64>j 


.3820 


.5325 .6492 




-.3780 


.3336 -.3113 




-.0292 


,0215 .0937 


a] 

is 


.1880 


.1285 .1009 


is 


.1582 


.1147 




-.0039 


.0576 




+.0019 


-.0418 


ih(jf} 

d^j 


+.3207 


.1296 -.3559 


mi 


-.0487 


.0467 .1813 


96y 


+ .2720 


.0829 -.1746 


S>03 


+ .0946 


-.1287 


%si:s 


-.0221 


+.0388 + 


«■ 


+ ■,00094 +.00068 + 





®3 


^3 


.3775 


.4164 


.4524 


.7399 


.8161 


.8867 


.3065 


-.3097 


-.3133 


.1768 


.2664 


.3619 


.0898 


.0652 


.0818 



.0953 


.0875 


.0835 


.1352 


.2216 


.3141 


.1308 


-.2501 


-.4043 



,4518 


-.4635 


-.4084 


.3356 


.5017 


.6786 


.1162 


+.0382 


+,2702 



.1454 


-.0390 


+.1542 


.0515 


+.0155 


-.0670 


.00008 


-.00066 


-.00158 
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The constant is found by applying Equation ( 77) at the inner 
channel boundary. 



^ s 3.7894-.14286[.1429+.0380+. 0009 ^-2(l. 3438-. 2873^ = 3.4616 

n n-1 

Pj= .5878(1.0443)^ = .7000 C2«(.7000) ^ =.9153 





^3 




*^3 


03 


®3 


^3 


n-1 


P 3 n 


.9153 


.9162 


.9169 


.9170 


.9163 


.9148 


P3 


.7000 


.7028 


.7049 


.7052 


.7031 


.6984 


^3 


1. 1909 


1. 1191 


1.0706 


1.0345 


1 . 0044 


.9773 


P- i 




1. 1404 


1.0883 


1.0526 


1.0258 


1.0033 


.9829 



2 



Av. 

n-1 


1. 1143 


1.0704 


1.0392 


1.0145 


.9931 


P ” 
^3 


.9157 


.9165 


.9169 


.9166 


.9155 


/ Ct p 71 


+.0380 


+.0227 


+.0026 


-. 0209 


-. 0490 



The axial velocity integral is iterated in steps across the channel as 
follows: 



<Ps 

. 02083[l. 8139(9. 1313-8. 846l)+3.92f. 0946 ^.078§= .09602 .09602 
. 01874[l. 8797(9. 0490-9. 0107)+3. 92-. 1287 K 048’^= .07331 .16933 
. 01716[l. 9278(8. 9897-9. 2287)f3.92-. 1454K005^= .05697 .22630 
.01593[l. 9610(8. 9174-9.3537)^3. 92-. 0390-.04s|= .04742 .27372 
.01495[l. 9811(8. 8593-9.4040) ^3,92+. 1542-. 117|= .04302 . 31674 



Average 
. 04801 

. 13267 

, 19781 

. 25001 

. 29523 



The mass flow equation is iterated to solve for as follows: 



I 






- 1 



I 



1 

I 
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.026865=. 017071(. 04801 H. 1820Cj)+. 016377(. 13267+1. 0766Cj) 
+. 015754(. 19781+. 9985Cj)+. 015167(. 25001 +. 9374Cj) 

+. 014597(. 29523+. 83770^) 



r- _ .012654 _ 1 C/.0 

^1 “ :7jmi7 - -1568 



Solution at station 3: (n = I. 33) 





^3 


‘^3 


^3 


<^3 


®3 


^3 




.7600 


.8661 


.9605 


1.0464 


1. 1258 


1. 2000 




. 1949 


.2717 


. 3312 


.3775 


.4164 


.4524 


^3 


-. 3780 


-.3336 


-.3113 


. 3065 


-.3097 


-.3133 


9^ 


-.0292 


. 0215 


.0937 


. 1768 


. 2664 


.3619 


P 3 


.7000 


.7028 


.7049 


.7052 


.7031 


.6984 


'^3 


3.4617 


3.6285 


3.7862 


3.9354 


4. 0740 


4. 2038 


O 3 


. 20220 


. 19367 


. 18615 


. 17921 


. 17275 


. 16518 




3.4876 


3.6503 


3.8062 


3.9582 


4. 1049 


4. 2456 


^3t 


.7185 


.7177 


.7181 


.7204 


.7219 


.7229 


^3t 


. 20600 


. 19659 


. 18863 


.18183 


. 17605 


. 16929 




-.0001 


-.0047 


-.0108 


-.0155 


-.0178 


-.0202 



Compared to the isentropic case (Examples Vla^b), the poly> 
tropic flow with n = 1. 33 through the same stator and rotor results in 
the following maximum percentage changes at station 3: 

+18.0®/o T^: •!. l°/o 

+20.8°/o ^3t:-9.3®/o 

Pj: - 8. 3*^/o Workout; +19. 5°/o 

Example IXa 

Solution for the stator of Fig. 5 (divergent curved channel 
walls) with a radial total temperature gradient and radial coordinates 
for equal areas. This is an isentropic, subsonic solution for the direct 
problem for the stator. The resulting solution at station 2 is used in 



8i 



Example IXb as the upstream conditions for a solution of the 
problem for a rotor. 



The equations used in Example Ilia are applicable. 
Conditions at station 1 : given 







^1 

.8000 


.8438 


^1 

.8854 


^1 

.9252 


®1 

.9633 






3.7894 


3.9679 


4. 1374 


4. 2995 


4.4548 


Pit 




1.00 


1.00 


1.00 


1.00 


1.00 


Qlt 




. 26389 


.25202 


.24170 


.23258 


. 22448 


Pi 




.9400 


.9400 


.9400 


.9400 


.9400 


r, 




3.7231 


3.8985 


4.0650 


4. 2243 


4. 3768 






. 25249 


.24113 


.23126 


. 22253 


.21478 


G, 




-.0340 


.0033 


.0403 


.0768 


. 1125 


4>, 




.6804 


.6970 


.7103 


.7215 


.7303 


F - 


2.75 


m = .029242 


dTit . 


0746 


Solution at station 2: 

“2 *>2 
^2 . 7850 . 8496 


.9098 


^2 

.9661 


®2 

1.0194 


<Pz 




.7008 


.6688 


.6452 


.(iZlS 


.6139 






1.9272 


1.8392 


1.7743 


1.7256 


1.6832 






-.0701 


.0241 


. 1050 


. 1765 


.2416 


P2 




.5460 


.5954 


.6322 


. 6601 


.6822 






3. 1379 


3.4217 


3.6292 


3.3183 


3.9938 


Oz 




. 17126 


. 17402 


. 17419 


. 17286 


. 17081 






0 


-.0010 


-.0026 


-.0049 


-.0083 


T - T 

~ *-/t 

Cj = . 5659 
Example IXb 


P2t 


'Pit’ ^ 


2t ' 
8412 





inverse 



0000 

6043 

00 

21719 

9400 

5237 

20781 

1473 

7365 

0 



0700 

6033 

6591 

3016 

6997 

1577 

16828 

0116 



Solution for the rotor of Fig. 5 (divergent curved channel 
walls) witli a radial temperature gradient and radial coordinates for 
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eq-iial areas. This is an isentropic, subsonic solution for the in- 
verse problem (H{^) given) for the rotor. One complete approxi- 
mation for the inverse problem for the rotor is given below to il- 
lustrate the iteration procedure. 

Equations (63), (65), (73), (75), (77) (with appropriate 
subscripts) emd the following two relations are used; 



(h - I ^ ^ (f 2^ 2\cLf/ cLds 



(80) 






(27). 



When iterated across the channel. Equation (80) becomes: 







(81). 




Conditions at station 2 are the same as after the stator, 
Example IXa above. The work output function H(^ ) is given in this 
example, and the blade shape (spouting angle) is determined. 



Data used in tlie iteration: 

diTzi 



m = . 029242 
(Hj-H.^P = 0 



dL4. 



2.8593 



cLH 

<^6 



= 0 



H(/^) = 1. 300 



(constant- 

given) 



I 
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^3 


'^3 


®^3 


^3 


®3 


^3 




.7600 


.8661 


.9605 


1 . 0464 


1.1258 


1.2000 




1.3000 


1.3000 


1. 3000 


1.3000 


1.3000 


1.3000 




1.5129 


1.5626 


1.6143 


1.6671 


1.7210 


1.7752 


^3 ^3 


-. 2129 


-;2626 


-.3143 


-.3671 


-.4210 


-.4752 




-.2801 


-.3032 


-. 3272 


-.3508 


-. 3739 


-.3960 



/ 





1.7889 


1.7317 


1.6766 


1.6238 


1.5736 


y cLtzt(4i-4^^ 
r-, ^ 


.86305 


.86305 


.86305 


.86305 


.86305 




-.02931 


-.02801 


-.02751 


-. 02681 


-.02631 




. 1061 


.0944 


.0859 


.0794 


.0742 


odf-c:) 
^2 2 


. 01489 


.01501 


.01496 


.01482 


.01462 


¥iS^rL) 


. 007820 


. 006658 


.005821 


.005185 


. 004682 




1 ' 
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Try Cj =.1039. 





as 63 


*3 


. 1539 . 2379 


A 5 


-. 2801 -.3032 


63 


-.0156 .0188 




.1047 .1038 




.0016 .0560 


903 


.3242 .7881 




-.0005 -.0441 


d4^3 


-2.9845 -2.1545 


0<p3(d-^) 

ipCJs 


+.4477 -.1704 


Mw), 


-. 0260 . 0409 


903 


+.4217 -.1295 




+.1461 -.1978 


-S^ 


-.0065 +.0592 



*^3 




®3 


^3 


3295 


.4159 


.5242 


.6891 


3272 


-.3508 


-.3739 


-.3960 


0932 


. 1948 


.3354 


. 5513 



. 1145 


. 1209 


. 1274 


. 1440 


.2651 


.4433 


1.1827 


1.7707 


2.9097 


-. 1703 


-.4694 


-1.2889 



.5785 


-1.0580 


-.4485 


+.4290 


. 4466 


-.4956 


-.2870 


+.3432 


1804 


.3698 


. 6316 


1.0336 


2662 


-. 1258 


+. 3446 


+1.3768 



-.I960 +. 1094 K8607 

+.0700 -.0642 -.5647 



7!d 
tlR ^ 



+.00076 +.00082 +.00008 



00214 



00809 



i 



I 



L 
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The constant C, is found by applying Equation (77) at the inner 
channel boundary. 

^ = 3. 7894 - . 14286 [. 0785 + . 0108 + . 0002 + 2. 60o] = 3. 4052 

■f 211 

P3 =.5460 (1.0682)»'-/ =.6878 C2 =(.6878)^ = 8986 



r-f 

P 3 ’' 


^3 

. 8986 


63 C 3 

.8994 .9002 


<^3 

.9003 


®3 

, 8981 


^3 

. 8900 


P 3 


. 6378 


. 6900 . 6921 


. 6924 


. 6865 


.6651 


) ‘ 


. 2597 


1.1589 1.0947 


1. 0489 


1. 0063 


.9506 




. 1793 


1.1111 1.0663 


1 . 0347 


1. 0045 


. 9645 


Av. 


1 . 


1452 1.0889 1. 


0507 1.0196 .9845 




(i-V) 


. 1014 


.1006 .0998 


. 0997 


. 1019 


. 1100 



Av. 



. 1010 



. 1002 



. 0998 



1008 



. 1059 



The axial velocity integral is iterated in steps across the 
channel as follows: 

^ <p3 

2. 0486 (. 08717 . 02931) + . 01550 = . 13403 

1. 8856 (. 08648 02801) - . 01867 = . 09158 

1. 7616 (. 08613 02751) - . 01684 = . 03642 

1. 6556 (.08700 - .02681) + .00869 = . 10834 

1. 5492 (. 09140 - . 02631) t. 06386 = . 16470 



029242 = .017052 (. 06701 + C^) + . 016344 (. 17982 + C^) 

+ .015718 (.26882 + Cj) + .015110 (.36620 + Cj) 

^ . 014393 (. 50272 + Cj) 

^ _ . 008166 _ , 

<=l"-07SFT7 





A.verage 


. 13403 


. 06701 


. 22561 


. 17982 


. 31203 


. 26882 


.42037 


. 36620 


. 58507 


. 50272 


solve for Cj as fo 




I 
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Solution at station 3; 





^3 


*^3 




<^3 


^3 


'3 




. 7600 


. 8661 


. 9605 


1. 0464 


1. 1258 


1. 2000 


<ps 


. 1039 


. 2379 


. 3295 


.4159 


. 5243 


. 6890 




2801 


-. 3032 


-. 3272 


-. 3508 


-. 3739 


-. 3960 


Os 


0156 


. 0188 


. 0932 


. 1948 


. 3354 


. 5512 


P 3 


. 6878 


. 6900 


. 6921 


. 6924 


. 6865 


.6651 


Ts 


3.4052 


3. 5688 


3. 7243 


3.8706 


4. 0010 


4. 0978 


^3 


. 20197 


. 19335 


. 18583 


. 17886 


. 17158 


. 16231 


G//J 4.619 


2. 366 


1.922 


1. 672 


1.434 


1. 167 


(O3 


12°12’ 


22°55' 


27°30' 


30°55' 


34°53' 


40°35' 




0 


+ . 0064 


+. 0096 


+. 0097 


+ . 0070 


+ . 0014 



Note: A supersonic solution is indicated for problem IXb by using 
average values in the mass flow and energy equations 
{ <f>^ = 3.780). However, the supersonic solution is im- 
possible in this case since the pressure goes to zero at the 
outer channel boundary before the solution is reached. 

Example Xa. 

Escamples Xa, b are the isentropic, subsonic solutions for the 
stator and rotor of Fig. 5. This is an example of the inverse problem 
in which H( ^3 ) and = 0 are specified and the blade shapes of 

both the stator and rotor are tlius determined. 

Equations (27), (63 ), (65), (73 ), (74), (75) (with appropriate 
subscripts) are applicable for the stator in addition to the following 



relation: 
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_l_ d.(^ z '^z) 

<P^ J 






( 82 ). 



Conditions at station 1 are assumed the same as for Example 



IXa. 



^2 


^2 


^2 


<^2 


®2 


^2 


1. 3000 


1. 3561 


1.4061 


1.4515 


1.4935 


1. 5328 


(given) 

i 1.6561 


1. 5962 


1. 5455 


1. 5024 


1.4651 


1.4325 



m = . 029242 



Solution at station 2: 



dr,t 



= 4. 0746 




Aa 

4>z 



^2 


^2 


^2 


*^2 


®2 


'2 


^2 . 7850 


.8496 


. 9098 


.9661 


1 . 0194 


1. 0700 


4>^ .6791 


.6327 


. 6004 


. 5760 


. 5532 


. 5311 


6561 


1. 5962 


1. 5455 


1. 5024 


1.4651 


1.4325 


9^ -. 0679 


. 0228 


. 0977 


. 1620 


. 2177 


. 2656 


P 2 .6367 


. 6780 


.6932 


.7308 


' . 7489 


.7635 


3.3310 


3. 5501 


3.7485 


3. 9310 


4. 1015 


4. 2627 


Q 2 .19116 


. 19093 


. 18604 


. 18590 


. 18279 


. 17912 


2. 439 


2. 523 


2. 574 


2. 608 


2.648 


2. 697 


22°18' 


21°36' 


21®15' 


20®59' 


20°41' 


20°21' 


0 


-. 0035 


-. 0052 


-. 0051 


-. 0038 


-. 0019 


Cj = .6791 C 2 


= . 8790 


^2t = ^If 


^2t "" ^'It. 


Z - T 
at It 



I 
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Example Xb. 



Solution for the rotor of Fig. 5 for the inverse problem in 
which H ( ^ ) and 0 are specified and the blade shape O 

is to be determined. 

Equations (27), (63 ), (65), (73), (75), ( 77) (with appropriate 
subscripts are applicable for the rotor in addition to the following 
relation : 




Conditions at station 2 are the same as after the stator, 



Example Xa. 



m = . 029242 


dTzi: _ ^ 


8593 


= 0 


Solution at station 3: 










^3 


63 


^3 


<^3 


®3 


'3 


. 7600 


. 8661 


.9605 


1 . 0464 


1. 1258 


1 . 2000 


05 . 2985 


. 3398 


.3627 


. 3830 


.4110 


.4602 


9s -.0048 


. 0269 


. 1026 


. 1794 


. 2630 


. 3682 


0 


0 


0 


0 


0 


0 


P 3 . 6876 


. 6873 


.6868 


. 6857 


. 6825 


.6758 


Q 3 20196 


. 19280 


. 18431 


. 17763 


. 17105 


. 16416 


T 3 3.4049 


3. 5639 


3.7388 


3.8602 


3. 9940 


4. 1165 


Clfj) 2 . 546 


2 . 549 


2.648 


2 . 732 


2.739 


2 . 608 


21®25' 


21°23‘ 


20°41' 


20°07' 


20 ^ 02 ’ 


20^57' 


t> 

0 


0001 


+.0011 


-. 0010 


+. 0001 


+. 0002 



Cj = . 2985 



C 2 = .8985 
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Eacample XI 

This example has been constructed to illustrate the increase 
in gross thrust obtainable by the use of a radial total temperature 
gradient. 

A hypothetical one-stage turbine consisting of a stator and a 
rotor has been assumed similar to that shown in Fig. 3. In order 
to give this example more physical meaning, rotor blades have been 
used which are tapered linearly in area in such a way that the ratio 
of hub to tip cross-sectional area is three. Using a rotor tip speed 
of 1100 feet per second, the stress due to centrifugal forces was 
calculated at various points along the blade. Using the high-tcmpera- 
ture data for ’'Vitailium”, (H.S. ?A) (NR-iO), given in References 16 
and 17, it was then estimated that the blade root could safely stand a 
temperature of 1450°F, while the blade tip temperature could be 
1850 ”f. 

Using this data, the thrust of this turbine is calculated for 
two conditions: (a) A constant radial total temperature of 1450°F, 
and (b) a linear temperature gradient of 400 F with the same root 
total temperature of 1450°F. The mass flow, boundary configuration, 
combustion chamber total pressure, and worl: extracted arc the same 
in both cases. 

Example XIa. 

Solution for the stator and rotor of Fig. 3 with no radial total 
temperature gradient and equal area coordinates. This is an isen- 
tropic, subsonic solution to the inverse problem (H ( ^ 3 ) given). 



I 
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For no radial total temperature gradient, H ( ^ ^ ~ 

= constant, and a cylindrical channel. Equation (82) reduces to 
= Cj . Equations (63), (64), (65) and (67) are also used in 
the solution for the stator . 

Conditions at station 1: (given) 



= 2.7071 


Fit = 1-00 


= .36940 


.4915 


Z = 2.6726 


Fj = .9562 


Qj = .35775 


o 

II 


= 880 ^ 
sec 


= 1910° 


R (const .) V/ = 


148.26 

sec 



<o/P=1100-^ p,. = 16,300 - 1 ^ K(/)= 1 .3964 (const.) 

^ £> C C X w *3 

/?o=1.25ft m =.031650 A_y=0 



Solution at station 2: 





^2 


^^2 


°2 


^*2 


®2 


/a 


.8000 


.0544 


.9055 


.9539 


1 .0000 


<P^ 


.7260 


.7260 


.7260 


.7260 


.7260 




1.7455 


1 .6342 


1.5421 


1 .4639 


1 .3964 


^2 


.4812 


.5236 


.5585 


.5879 


.6125 


^2 


.21909 


.23268 


,24369 


.25277 


.26027 




2.1965 


2.2501 


2.2921 


2.3259 


2.3533 


SS ~ 


^ 2.404 

<P2 


2 .251 


2.124 


2.016 


1 .923 


/^e 


22°34' 


23°59' 


25°13' 


26°25‘ 


27°30' 


^‘2 


1 .071 


1 .005 


.951 


.904 


.873 




0 


+.0002 


0 


- .0003 


-.0001 




Cl = 


.7620 


II 

Si. 


^2t = "^It* 


^2t " ^It 




^2 = 


.8114 









8114 




I 
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For the rotor with no radial total temperature gradient, 

K ^^ 3 ^ ~ constant, and =0, the solution at station 3 

is easily found to be: 





- . 7806 


C3 = . 22525 










= 0 


= 2. 3081 








r> 

' 3 


= . 5003 


= . 24798 








^3 


= 2.2211 


P3t = -5724 












^3 ^3 


^3 


^3 


®3 




,«i I* 

v-e 

II 


1. 025 1.095 


1. 160 


1. 222 


1.281 






44^20' 42°25' 


40°46' 


39°17' 


37°56 



Assuming isentropic expansion to atmospheric pressure through 
an ideal nozzle, the total gross thrust is found to be 1 1 , 966 pounds . 
Example Xlb . 

Solution for the same stator and rotor as Example XIa except 
with a linear radial total temperature gradient of 400°F. The mass 
flow, combustion chamber total pressure, and work extracted are the 
same as in the above example. 

Equations (27), (63), (64), (65), (67), and ( 82 ) are applicable. 
Conditions at station 1: (given) 





ai 




^1 


^1 


®1 




. 8000 


. 8544 


.9055 


.9539 


1. 0000 


r,t 


2 . 7071 


2.8614 


3. 0062 


3. 1434 


3. 2741 


^it 


1. 00 


1. 00 


1. 00 


1. 00 


1.00 


^it 


. 36940 


. 34948 


. 33265 


. 31813 


. 30543 


^1 


. 9512 


.9512 


.9512 


. 9512 


.9512 


ry 


2.6687 


2.8208 


2.9635 


3. 0988 


3 . 2276 


<>i 


. 35643 


. 33721 


. >2097 


. 30696 


. 29471 


4>, 


. 5185 


. 5331 


. 5467 


. 5587 


. 5705 




0 


0 


0 


0 


0 
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m = . 031650 
Solution at s 


tation 2 : 
^2 


^^'^= 2.8350 H(4^) 

^2 


= 1. 3964 (const .), A^= 
^2 ®2 


4>a 


. 5844 


. 6975 


.7955 


. 8819 


.9589 


Aa 


1. 7455 


1. 6342 


1. 5421 


1. 4639 


1.3964 


?2 


. 5018 


. 5443 


. 5774 


. 6036 


. 6249 


U 


2. 2231 


2.4050 


2. 5697 


2. 7212 


2. 8625 


^2 


. 22574 . 22632 


. 22474 


. 22183 


. 21829 


^ ' 4>a 


\ 

2.987 


2. 342 


1.938 


1.659 


1 . 467 




18^30' 


23°07' 


27® 18' 


31°09' 


34° 19' 


lvl2 


1 . 044 


.974 


.920 


.878 


. 850 




0001 


+. 0054 


f. 0064 


+. 0049 


t. 0017 


= 

^2 = 


. 5844 
. 8212 




• ^2t 


~ ^It’ ^2t 


“ ^It 



For the rotor with h. (f^) = constant and - 0, 

Equation (64) reduces to = constaiit. Equations (63), (65), ( 71), 
and the follov/ing integral are applicable; 



i-i) 7.^, 



(84). 
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Solution at station 3; 







*^3 


*^3 


®3 


.6002 


.7426 


.8514 


.9397 


1 .0141 


)\ j 0 


0 


0 


0 


0 


^3 .5289 


.5289 


,5289 


,5289 


.5289 


Tj 2.2566 


2.3853 


2,5060 


2.6203 


2.7293 


Q- .23435 

3 


.22171 


.21103 


.20182 


.19376 


GO 

O 

• 


2 .4624 


2.6072 


2 . 7444 


2.8751 


F3, .5724 


.5912 


.6075 


.6220 


,6345 


.24798 


.24005 


.23300 


.22658 


,22067 


G = ^ 1 . 333 


1.151 


1 ,064 


1 .015 


.986 


Ps 36°52 ' 


40°58' 


43°13' 


44°35' 


45°22' 


0 


- .0017 


-.0023 


-.0020 


-.0011 




Cj = .6002 




= .8336 





Assuryiing isentropic expansion to atmospheric pressure 
through an ideal nozzle, and summing the thrust by stream tubes, 
the total gross thrust is found to be 12,883.8 pounds. This is 
an increase of 917,8 pounds, or 7.7 percent, over than found in 
Example XIa. 
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^3, ^3 


a 


.8000 


b 


.8544 


c 


.9055 


d 


.9539 


e 


1 .0000 



0.90 



0.80 



0.70 



0.60 



P 

d| 






1 

( 


4 






I ^ — — 1 
d3 


C| 


Stol 


or 


1 


TV i 

R( 

C2 J 


Dtor 




C3 




b| 








b2 ^ 


lamlines 




^3 


Q| 


_ 


> 




4 






03 


XWN 




/ 

/ 

/ 

/ 


. \\\^ 


k \\\ 


Station 1 




^ station 2 




^station 3 


> 


s 


> 




> 

► 



0.10 0.20 

Axial Coordinate 



0.30 0.40 0.50 



FIG.3- CHANNEL CONFIGURATION AND COORDINATES FOR 
EXAMPLES I, H, AND 21 
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^2 


^3 


' ^ ^ '1,2,3 


a 


.800 


.790 


.780 


- .050 


b 


.850 


.850 


.850 


0 


c 


.900 


.910 


.920 


+ .050 


d 


.950 


.970 


.990 


+ .100 


e 


1 .000 


1 .030 


1 .060 


+ .150 




Axiol Coordinate 



FIG. 4 — CHANNEL CONFIGURATION AND COORDINATES FOR 



EXAMPLE nr 




i 




Radiol Coordinate 
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Axial Coordinate 

FIG. 5 - CHANNEL CONFIGURATION AND COORDINATES FOR 



EXAMPLES 12-3: 



I 
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Exomples 13E , X 





s , 


^ 2 


^ 3 


/Jl) 

Us/| 








Q 


.800 


.785 


.760 


-.05 


-.10 


-.15 


-.25 


b 


.840 


.842 


.848 


0 


.02 


,04 


.10 


c 


.880 


.899 


.936 


.05 


.14 


.23 


.45 


d 


.920 


.955 


1 .024 


.10 


.26 


.42 


.80 


e 


.960 


1 .013 


1.112 


.15 


.38 


.61 


1.15 


f 


1 .000 


1 .070 


1 .200 


.20 


.50 


.80 


1.50 



Examples 3ZI -X 









4s 


(tI), 


(lil 


Oil 


(S4 




84 


a 


.8000 


.7850 


.7600 


- .0500 


-.1000 


- .1500 


- .2500 


-.2500 


- .2500 


b 


.8438 


.8496 


.8661 


.0047 


.0350 


'.0791 


.1332 


.1467 


.1720 


c 


.8854 


.9098 


.9605 


,0567 


.1627 


.2829 


.4972 


.5163 


.5474 


d 


.92 52 


.9561 


1 .0464 


.1055 


.2813 


.4684 


.8455 


.8620 


.8891 


e 


.9633 


1 .0194 


1 .1258 


.1541 


.3935 


.6398 


1 .1789 


1.1893 


1 .2049 


f 


1 .0000 


1 .0700 


1 .2000 


.2000 


.5000 


.8000 


1 .5000 


1.5000 


1 .5000 



TABLE X - STREAMSURFACE COORDINATES, SLOPES AND CURVATURES 
USED FOR EXAMPLES JZ -X AND FIGURE 5 
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ROTOR FROM SUBSONIC STATOR 
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I 
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Radiol Coordinote 
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I 

1 

I 

1 
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FIG. 8 -TYPICAL TURBINE VELOCITY DIAGRAM SHOWING SIGN 
CONVENTION USED IN ALL EXAMPLES 
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^3 3 



Stator: § g) = ? 3) = Rotor: Cot /53 ( ^ 3) =0(13)= 



FIG . 9 -TYPICAL COMPRESSOR VELOCITY DIAGRAM SHOWING 
SIGN CONVENTION 
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